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Individual Elemental Abundances in Elliptical Galaxies 

Abstract 
by Jedidiah Lee Serven, Ph.D. 
Washington State University 
August 2010 



Using synthetic spectra to gauge the observational consequences of altering the abun- 
dance of individual elements, I determine the observability of new Lick IDS style indices 
designed to target individual elements. Then using these new indices and single stellar pop- 
O ' ulation models, I investigate a new method to determine Balmer series emission in a Sloan 
Digital Sky Surveys grand average of quiescent galaxies. I also investigate the effects of an 
■ old metal-poor stellar population on the near ultra violet spectrum through the use of these 
I ■ new indices and find that the presence of a small old metal-poor population accounts for 
^ ■ discrepancies observed between index trends in the near UV and optical spectral regimes. 
^ ■ Index trends for 74 indices and three data sets are presented and discussed. Finally, I de- 
termine the near nuclear line-strength gradients of 18 red sequence elliptical Virgo cluster 
galaxies for 74 indices. 
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CHAPTER ONE 
INTRODUCTION 

The question of how galaxies form is a far from simple matter. Recent results seem to 
make it look simple, in the young universe (z ~ 6) to the middle aged universe ( z ~ 1.5) 



the st ar formation rate for galaxies is approximately constant and then falls off ( iHartwick 
( I2OO4J ) and references therein). However, a closer look at the properties of galaxies such as 
morphological type and stellar content reveals a complex evolutionary picture. Just some of 
the parameters that an evolutionary picture must fit are (1) spirals must have formed from 
an accretion disk of gas. Ellipticals, on the other hand, may have formed from any number 
of scenarios, including mergers of various numbers and types of progenitor galaxies. (2) 
There is a scatt er in the mean ages of elliptical galaxies, which is less for larger galaxies 
( 1Wortheyll2007l ). And, (3) there is a scatter in the elemental abundance ratios from galaxy 
to galaxy and a systematic enhancement in large ellipticals for light elements such as Mg, N 
and Na. The most likely cause of this element enhancement is a change in the ratio of Ty pe 
I to Type II supernovae that are responsible for the chemical enrichment (jWortheylll998l ). 



The chemical makeup of these stellar systems can be studied using an integrated light 
spectrum. These integrated-light spectra formed from the added contributions of all the 
stellar light captured by an observer in the absence of complications such as emission and 
dust screening are indicators of the ages and elemental abundances of the stars in that 
system. Stellar systems such as elliptical galaxies, SO galaxies, star clusters, and the bulges 
of spiral galax ie s tend to be relatively free of gas and dust compared to disk galaxies 
( iKnappI ( I1999I ): lOosterloo et al.l ( 120021 )). These types of stellar systems are likely to contain 
small amounts of young stars, nebular gas, or dust making their spectra useful for 
comparison to stellar population models. For this thesis, a stellar population is defined as a 
single stellar population ( SSP ) of one age and composition, but all stellar masses. "Stellar 
populations" for a galaxy would be the weighted, collection of N-dimensional voxels, where 
the dimensions are age, average metallicity (<Z>), and the various individual metallicities 
<Z,>. 

Current models agree that age and abunda nce contro l the c o lor and spectral fea t ure 
strength of simp le integrated-light sv stems (IWortheyl (Il994j) : iBruzual &: CharlotI (|2003[ ): 
MarastonI (120051 )). The "3/2 rule" of IWorthey et al.l (jl994b[ ) states that if an age change is 
opposed by an abundance change such that d log (age) = -{3/2)d log (abundances), then 
very little change in colors, s urface brightne s s fluctu ation magnitudes, and absorption 
feature strengths will occur. IWorthey et al.l ( 1l994bl ) pointed out that this could be used to 
discover non-lockstep abundance ratio variations for individual elements, without having to 
know the absolute age or metallicity to very high precision if features could be found in the 
spectrum that were sensitive to such non-lockstep behavior. 
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Worthey et al.l (119921 ) presented models and observations that indicated that individual 
elemental abundances for individual galaxies could be found. This was shown by first 
plotting Mg and Fe absorption features for both models and observed elliptical galaxies. 
The models were for different ages and metallicity, and fell on a one-dimensional, highly 
degenerate trajectory where any reasonable change in age or metallicity simply moved an 
object back and forth along this line. The galaxies on the other hand showed considerable 
scatter about these model lines and showed that Mg features were stronger for larger 
galaxies. This indicated that there wer e nonsola r abun d ance ratios. The presenc e of this 
effect has been confi r nied many times (IWortheyl ( 19981): iTantalo fc Chiosil (|2004[ ) ; 
Thomas et aP feoosh : IProctor fc SansonJ J2002h : Maraston et all tooi )). 



Trager et al.l ( 120001 ) introduced the idea of using synthetic spectra to simulate the effects of 
elemental abundance ratio changes in integrated-light models. The technique involves 
generating many models of varying age, overall metallicity and individual abundance ratios 
then comparing those models with an observed galaxy to find a best fit model. This is 
commonly done by comparing measurements of spectra indices and choosing the model 
that best fits the observed galaxy's index measurements. This best fit model is then used 
to determine the age, metallicity and various elemental abundances. This technique is in 
common use today by those studying formation histories of elliptical galaxies. 



Serven et al.l ( 120051 ) employed this technique of using synthetic spectra to simulate the 
effects of elemental abundance ratio changes to expand the possible use of model index 
fitting to as many as 23 elements by using a crude model galaxy consisting of stellar 
models for a turnoff dwarf and a giant then combing them using equal contributions from 
the turnoff dwarf and giant at 5000 A to simulate a 5 Gyr stellar population. Then, by 
comparing this model with itself except with the abundance of a given element doubled, for 
example, the abundance of Cr would be increased by 0.3 dex, new indices were defined that 
were designed to be sensitive to this element . 

The rest of this thesis explores the application of this synthetic spectra technique to 
determine elliptical galaxies characteristics such as ages, elemental abundances and 
elemental distributions. 



Chapter 2 presents the IServen et al.l (120051 ) expansion of index definitions to cover 23 
elements and the results of that work. 

Chapter 3 illustrates the use of index measurements of t he hydrogen Balmer series and a 
new technique to determine emission corrections for the I Graves et al.l (120071 ) grand 
averages of elliptical galaxies form the Sloan Digital Sky Survey (SDSS). 



Chapter 4 explores the use of these targeted indices to investigate discrepancies between 



6 



index trends in indices in the near-ultraviolet and optical spectral regimes. 

Chapter 5 presents index trends for 74 indices and three data sets. 

Chapter 6 presents the near nuclear line-strength gradients for 74 indices and 18 Virgo 
Cluster galaxies. 
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CHAPTER TWO 

The Observability of Abundance Ratio Effects in Dynamically Hot Stellar Systems 



Introduction 



This chapter represents an update of IServen et al.l ( l2005l ). The purpose of that paper was 
to investigate the possibihty of determining the abundances of C, N, O, Na, Mg, Al, Si, S, 
K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Ba, and Eu in real galaxies by using 
synthetic spectra to define absorption feature indices sensitive to the individual abundances 
of these 23 elements. The methods used for that investigation are discussed in §2. Then, §3 
discusses the updates used in this work. 



2. Analysis 



In IServen et al.l (120051 ) to find places in the spectrum that show a noticeable response to 
the abundance change of an individual element, a "model" galaxy was constructed by 
combining two stellar synthetic spectra. One was the synthetic spectrum of a giant {Teff = 
4250, log g = 1.90) and the other main sequence turnoff dwarf {Teff = 6200, log g = 4.10). 
These two spectra were then combined with equal contributions at 5000 A. The composite 
spectrum was then broadened to a simulated line-of-sight velocity field that is a Gaussian 
distribution with a width (a) of 200 km s~^ to simulate an elliptical galaxy a little more 
massive than the Milky Way. 

First this was done for scaled solar metallicity, then the same model was computed with 
the exception that an individual elemental abundance was increased by 0.3 dex. This was 
done for 22 of the 23 elements. The exception was C, which was increased by 0.15 dex due 
to the fact that for a 0.3 dex increase in C the C/0 ratio approaches 1 and thus drastically 
changes the structure of the star toward that of a carbon star. Finally to locate places of 
significant response the ratio of an enhanced galaxy and the scaled solar galaxy was taken. 
This produced a spectrum where changes due to the element under study would show up 
as deviations from one in the ratio. 

After finding these areas of spectral response. Lick IDS (Image Dissector Scanner) style 
indices were defined consisting of a central passband that spans the absorption feature 
itself and two continuum passbands, one on the blue side of the central band and one to 
the red. From these three passbands one can find a pseudo-equivalent width for any feature 
of interest. To do this one first constructs a pseudo-continuum by "drawing" a line between 



9 



the wavelength midpoints and average flux values of the blue and red passbands. Then by 
integrating the difference between the pseudo-continuum and the galaxy flux over the 
central or index passband one gets an equivalent width characterizing the feature of 
interest (See Figure [3TT1) . 



Co Enhancement and Index Definition 




7S00 

Wavelength in Angstroms 



Fig. 2.1. — An example of a Lick IDS style index is plotted. The blue continuum passband 
(blue), the red continuum passband (red), the central passband represented by the base of 
the trapezoid (black) and a model spectrum (black) are plotted. The pseudo-continuum line 
connecting the centers of the blue and red continuum passbands, which forms the top of the 
trapezoid (black) is also plotted. The area within the trapezoid and bound by the spectrum 
represents the difference between the pseudo-continuum and the galaxy flux. Integrating 
this difference gives the equivalent width characterizing the feature of interest. 



Using this method IServen et al.l (120051 ) defined 53 new indices for these 23 elements. The 
feasibility of measuring these elements using these indices was checked by determining 
index values and errors of these models and assuming a photon error such that S/N = 100 
around 5000 A and proportional to F^j^ . They determined that of the 23 elements 18 



looked like they could be measured and 5 ( S, K, Cu, Zn, and Eu 
determine. 



looked difficult to 



Discussion 



Serven et al.l ( 120051 ) then went on to determine the spectral response of these 53 new 
indices and 37 previously defined indices to these 23 elements shedding light on the 
infiuence any particular element may have on any given index. They summed up this work 
in two large tables giving the index response in angstroms of equivalent width. Since the 
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rest of this paper will often be referring to these indices and their sensitivities those two 
tables have been recreated in the form of 6 tables located in Appendix A. In these tables 
the 1st column of the table gives the indices's names, the 2nd gives the index measurements 
(/o) in angstroms of equivalent width, and the 3rd giving the error (cr) associated with S/N 
= 500 at 5000 A. The next 23 columns split between tables are the changes (enhanced 
minus unenhanced) in the index due to the elemental enhancement 0.3 dex (0.15 dex for C) 
of the element labeling the column in units of angstroms. The 24th column is the change 
brought about in the index by increasing all the alpha elements by a factor of 2. 

One major difference between these two works is that some of the indices that are not used 
in this work have been left out. The index definitions t hat are used include th e original 25 
l ick in dices ( jWortheyi Il994l : IWorthev fc Ottavianilll997l) , 1 Ha index defi ned in ICohen et al. 
( 1998 ). the 53 indices defined in IServen et al. f l2005h . 2 indices defined in iDavidge &: ClarkI 
(I1994J ) and 3 defined in chapter 3. 



The other major difference is that the numbers tabulated here come from measur ing these 



(1994 


) & 


Traeer et al. 


(1998) 



optical ( iLee et al.ll2009[ ) in order to investigate the effects of changing the detailed 
elemental composition on an integrated spectrum. These models where then used to create 
synthetic spectra at a variety of ages and metallic ities for si r igle-b urst stellar populations. 
The underlying isochrones for this paper were the IWortheyl (Il994| ) isochrones, because they 
allow us "manual" HB (Horizontal Branch) morphology control. However, there are certain 
caveats to using these isochrones. Specifically the models are a bit crude by today's 
standards and the ages are about 2 Gyr too old, so that 17 Gyr should really be 
interpreted as 15 Gyr. Other isochrone sets were used to check the results. 

For this thesis new index fitting functions were generat ed. The data sources i nclude a 
variant of the original Lick collection of stellar spectra ( iWorthey et al.lll994bl ) in which the 
wavelength scale of each observation has been refined via cross-correlation, as well as the 
MILES (Medium-resolution Isaac Newt on Telescope library of empirical spectra) spectral 
library ISanchez-Blazguez et al.l (l2006b) wi th some zero point corrections, and the Coude 
Feed Library (CFL) of IValdes et al.l (120041 ). The CFL was used as the fiducial set, in the 
sense that any zero point shifts between libraries were corrected to agree with the CFL 
case. The MILES and CFL spectra were smoothed to a common Gaussian smoothing 
corresponding to 200 km s~^. The rectified-Lick spectra were measured and then a linear 
transformation was applied to put it on the fiducial system. 

Multivariate polynomial fitting was done in five overlapping temperature swaths as a 
function of Gg// = 5040/Tejj, log g, and [Fe/H]. The fits were combined into a lookup 
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table for final use. As in IWortheyi ( 119941 ). an index was looked up for each "star" in the 
isochrone and decomposed into "index" and "continuum" fiuxes, which added, then 
re-formed into an index representing the final, integrated value after the summation. This 
gives us empirical synthetic spectra when variations in chemical composition are needed. 
The grid of synthetic spectra is complete enough to predict nearly arbitrary composition. 
For the remainder of this paper whenever the term model is used it is to these models that 
we are referring as they are the backbone of methods used in the following chapters. 
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CHAPTER THREE 

Emission Corrections for Hydrogen Features of the Graves et al. 2007 Sloan Digital Sky 
Survey Averages of Early Type, Non-LINER Galaxies 



Introduction 



One of the more important factors in determining the age of a stellar population such as an 
elliptical galaxy is being able to measure the hydrogen Balmer series (IWortheylll994 : 
O ' Connelll 1 1 9 761 ) . The first two Balmer lines Ha and H/3 are primarily used to determine 
the age of a stellar population because of their rela tive insensitivity to changes in 
metallicity ( IServen et al.l ( l2005l ): iKorn et al.l (120051 )) making a measurement of these two 
lines a more reliable estimator of age t hen the higher order Balmer lines H 7 and B.6, which 
are sensitive to changes in metallicity ( IServen et al.ll2005l : iKorn et al.ll2005l ) making a 
knowledge of the relative abundances of elements in a given spectra necessary in order to 
get a consistent agreement of ages from all the Balmer lines. 

Another important factor in determining the age of a galaxy is hydrogen emission, which 
can give the impression of a much older galaxy by filling in and weakening the diagnostic 
Balmer absorption features. In the past, an emission correction is used that is based on the 
on or OIII emission lines (jSchiavon et al.ll2006l : lGonzalezlll993[ ) in which the difference in 
the measured equivalent width (EW) of a hydrogen line such as Ha or H/3 due to emission 
is some fraction of the EW of Oil or OIII. However, the correlation between hydrogen and 
oxygen emission strengths is weak, and indeed there is no astrophysical reason why the two 
should be tightly correlated. On the other hand, being able to use hydrogen recombination 
l ines, which have ratios that are nearly fixed with respect to themselves 
( jOsterbrock Sz Ferlandll2006l ) should give a more reliable result. 



Toward that end, we compare 13 high-quality spectra of Virgo cluster ea rly-type galax ies 



(Serve n & Worthey, in preparation) with a sample of SDSS spectra from iGraves et al. 
j2007h . 



The Virgo Cluster elliptical galaxies consist of 13 red galaxies with colors 
0.75 <B-V< 0.97. They also have a high S/N from S/N=150 to S/N=450 per pixel, are 
well fluxed using standard star flux calibrations, and have tight control over instrument 
resolution and velocity dispersion. Long slit spectra were obtained using the T2KB chip on 
the Cassegrain spectrograph on the Kitt Peak Mayall 4 meter telescope (2006 Jan 31 - Feb 
5). A resolution of 1.4 A per pixel was selected to adequately sample the velocity 
dispersions of the program galaxies, which range from 50 to 350 km s~^. The wavelength 
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range covers from 3200 A to 7500 A in two wavelength swaths. 



The iGraves et al.l (120071 ) sample was constructed from spectra of 22,501 SDSS galaxies that 
fall into the redshift range 0.06< z <0.08 and whose colors ni eet the criteria (°-^)(g-r) > 
-0.025(°-^M^ - 51og h ) + 0.42, with h = 0.70 fiYan et al.l[2006[ l placing them firmly within 
the red sequence. These spectra are then further divided into those with Ha and Oil 
emission lines and those without. Those with emission lines and high [Oil] /Ha ratios were 
said to be "LINER-like" and those without "quiescent" where LINER stands for Low 
Ionization Nucle ar Emission-line Region. The high [OH] /Ha ratios are defined in 
Yan et al.l tood l as those with (EW([OII]) > 5EW(Ha) - 7). 



Graves et al.l (120071 ) found that the [Oil] equivalent width distribution of their quiescent 
sample has a standard deviation of a = 1.56 A. In part to reduce the impact of outliers 
which may contain emission just below the detection threshold, they produced a random 
sub-sample of 2000 quiescent galaxies selected to conform to a Gaussian distribution 
centered on an [Oil] equivalent width of and truncated at ±2cr[0II]. They then divided 
these 2000 galaxies into the following velocity bins: 70< a <120 km s~^, 120< a <145 km 
s-^ 145< a <165 km s-\ 165< a <190 km s-\ 190< a <220 km s-\ and 220< a <300 
km s~^. For uniformity in comparison, each spectrum in every bin was smoothed to an 
equivalent velocity dispersion of a = 300 km s~^, the highest velocity dispersion in the 
sample. Finally, they coadded all of the individual spectra in each bin to produce six 
composite spectra corresponding to quiescent galaxies with different original dispersions. 
Along with these spectra the S/N at each resolution element was computed to produce an 
error spectrum for each of the six composite spectra. Factors that may contribute to the 
error spectrum estimate are age, individual abundances, and emission signal. The error 
spectrum is likely to be dominated almost entirely by measurement uncertainty given the 
poor S/N of the individual SDSS spectra. 



We adopt the six composite spectra from lGraves et al.l (120071 ) for analysis and comparison 
with the Virgo spectra. The SDSS data points in the figures of this paper represent indices 
measured from those composite spectra. 

The central aim of this chapter is to measure the strength of any residual hydrogen 
emission in these quiescent spectra. Noting the apparent residual Ha emission of the SDSS 
spectra as seen in the first panel of Figure 3.1, we construct simple formulae for the 
correction of the relative intensities of the rest of the Balmer series for the Sloan Digital 
Sky Survey (SDSS) as a function of Mg b (a Lick Style Index measurement for Mg) that 
includes continuum slope effects and the intrinsic decrement values for the Balmer lines. 
The corrections should find future use with integrated-light models to better predict stellar 
populations parameters, especially the mean metallicities and ages of galaxies. 
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Our analysis methods and results are set out in the following section with a discussion and 
summary in § 3. 



5. Analysis and Results 



For this work, the Ha index is as defined by lCohen et al.l ( 119981 ). The first panel of Figure 



3.1 shows that there is a discrepancy in the Ha measurements between the SDSS and 
Virgo spectra. Plausibly, this is either due to hydrogen emission in the SDSS spectra or to 
nonspectrophotometric wavelength-dependent fiux calibration errors that would propagate 
into Lick index measurements depending on how the various passbands lay across the 
various spectral "lumps." We infer that the bulk of the discrepancy is most likely due to 
hydrogen emission in the SDSS spectra because the continuum-shape differences between 
the data sets are too small to account for much of the discrepancy, as we now show. 

To determine if differences in the continuum of the two data sets could explain this 
discrepancy, the ratio of similar spectra (similar as regards velocity dispersion) from the 
two data sets was computed and a continuum fit was found for this ratio using the 
"continuum" routine in IRAF 0. Lick style indices were measured from the fit as one would 
measure spectra. The difference in the index thus discovered was found to be small: about 
0.05 A, insufficient to explain the SDSS-Virgo discrepancy. It should be noted that the 
effects of the continuum differences could be more severe if differences in the continuum 
shapes exist that are of the same wavelength span as the features of interest. Most would 
likely lay with the SDSS spectra, since our fluxing of the Virgo data was careful, but 
unfortunately we do not have the data to check the signiflcance of these differences. 

We characterize d the SDSS-Virgo Ha - Mg b trend by best flt lines calculated using the 



routine fitexy.f ( iPress et al.lll992l ). This is a fortran program for finding the best fit line for 
data with errors in both the x and y coordinate. This routine minimizes the distance of 
each point from the line while taking into account weighting by the precision of the 
individual measurements in both the x and y coordinates. The choice to characterize the 
trends in Ha as functions of Mg b was made due to the fact that it is easy to model 
galaxies of various Mg abundance and easy to measure the strong Mg b feature. Also, since 
there exists a tight correlation between Mg b and a (velocity dispersion) these conclusions 
can be implied for galaxies of various a. Those line fits and the root mean square (RMS) of 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under Cooperative agreement with the National 
Science Foundation. 
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the distances of the points from their fit hnes are hsted in Table 3.1. Also shown in Table 
3.1 is the form of the correction term (jr^) determined from the following derivation. 



Table 3.1 



Data set 


6 

Line Fit in A 




RMS of Fit in A 


Virgo 


3.0132 - 0.3768*Mg 


b 


0.066 


SDSS 


1.3476 - 0.0532*Mg 


b 


0.020 


Correction Term 


1.6656 - 0.3236*Mg b -- 




0.069 



Table 3.1: Shown in this table are the best fit Ha vs. Mg h lines for the Virgo and Sloan 
data sets and their associated RMS values (lines 1 and 2). Line 3 is the difference between 
these two best fit lines, which represents the Ha correction term along with its RMS value. 



Using the linear correction term for the hydrogen emission in Ha a correction for the 
subsequent Balmer lines was constructed under the assumption that the entire shift is due 
to Balmer emission fill-in. To determine the form of the correction term we start with the 
definition of equivalent width (EW; see Eq. 1). In Eq. 1, Ai and A2 are defined as the blue 
and red wavelength bounds of the index passband, Fi{\) is t he flux at each wavelen gth 



across the passband, and Fc is the pseudo continuum flux; see IWorthey et al.l (119941 ). 



EW = - ^)d\ = AA(1 - ^) (1) 

t/ Al c c 

With more generality, and including a term for the flux due to an emission feature, 

EW = Td - fM±iM)dA = AA(1 - f - S) (2) 

where Fj is the flux of the emission feature and Fg is the flux of the stellar light and 

A A = A2 — Al. If we call the emission line's power j then the average emission line flux is 

defined by 

/•A2 

j= Fj{X)dX = AXTj (3) 

Thus, the average stellar flux inside the continuum band is Fg and the correction term for 
the equivalent width is AXFj/Fc = j/Fc. This gives us the equation for the equivalent 
width of an index with emission as 
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EW = rV - ^^^Il)d\ = AA - AA§ - ^ 
J\i Fc Fc Fc 

To extend to Balmer lines other than Ha, we exploit the fact that the decrements jp/j, 
j-y/ja, and js/ja are known, and relatively constant. 

For example, if ja is known, then extending from an Ha index to an H/3 index is 
accomplished by adding the correction term 



Fc,p 



X 



:,/3 ia 



(5) 



Having determined the Ha correction from best line fits (see Table 3.1) all that is left to do 
is determine the conversion factors and ^. 

Models (chapter 2 section 3) were used to determine the continuum level conversion 
factors(^^). The first step was to take the ratio of the continuum levels near Ha with the 
continuum levels near the other Balmer lines as measured from these models. The next 
step is to plot these continuum ratios against Mg h. The last step is to fit a least squares 
line to this data, giving a prescription for the change of the relative continuum levels as a 
function of Mg h. 

The last piece of the puzzle is that the native relative intensities of the Balmer lines need 
to be account ed for. This is taken into acc ount by the relative Balmer line intensities as 
calculated by lOsterbrock fc Ferlandl ( 120061 ). We adopt case B, 10000 K conditions, because 
they are near the middle of the range for star formation regions {ja/ jp = 2.85 ), and they 
are no t too drastically different than LINER type spectra (ja/j/? = 3.27) from lGraves et al. 
( 120071 ). Another reason to use these theoretical decrements is that observed decrements 
may suffer from asymmetric systematic errors due to local du st. The final correction 
formulae were apphed to H/3 defined in IWorthey et al.l ( 1l994l ). H7A, H7J7, H(5a, and ][i5F as 
defined in IWorthey fc Ottavianil (119971 ) and are found in Table 3.2. 
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Table 3.2 



Balmer 


Correction 


H« Correction 




Continuum Correction 


Decrement 


Correction term 


Index 


Term in A 










(^) 


Uncertainty in A 


H/3 


_ 


(1.666 - 0.324(Mg b)) 


X 


(0.838 - 


h 0.076(Mg b)) X 


0.351 


0.076 


H7A 




(1.666 - 0.324(Mg b)) 


X 


(0.847 - 


h0.149(Mg b)) X 


0.165 


0.093 


H7f 


j-f^ _ 


(1.666 - 0.324(Mg b)) 


X 


(0.847 - 


|-0.149(Mg b)) X 


0.165 


0.093 


USa 


if — 


(1.666 - 0.324(Mg 6)) 


X 


(0.649 - 


h 0.293(Mg 6)) X 


0.091 


0.140 


KSf 


is _ 


(1.666 ~ 0.324(Mg b)) 


X 


(0.649 - 


h0.293(Mg b)) X 


0.091 


0.140 



Table 3.2: The most reliable results are obtained for 3.0 A < Mg 6 < 4.3 A. This Mg h 
range covers most elliptical galaxies. The correction term uncertainty is calculated by the 
propagation of errors for each of the three terms in the correction formulae. The first error 
term is that of the Ha correction term found in Table 3.1. The second is the RMS of the 
continuum corrections term fit calculated in the same manner as that of the Ha correction 
and the third is due to scatter in the possible theoretical values of the decrement. 



Figure 3.1 shows the measurements of Virgo and SDSS galaxy averages along with the 
corrected SDSS data (light blue line). In all the graphs the Virgo data is in red while the 
SDSS data is in green. Also in Figure 3.1 model grids are plotted in blue and pink. The 
blue corresponds to models of various ages from 1.5 to 17 Gyrs and the pink corresponds to 
models of various metallicities from -2.00 to 0.50. For the H/3 correction the line fit is a 
little low for smaller galaxies with weaker Mg 6, but not bad for the larger galaxies (See 
Figure EH). 



18 




Fig. 3.1. — In panel 1 Ha, H/3, and H7A against Mg b indices for two data sets and models. 
In panel 2 H7^, B.Sa, and B.6f against Mg b indices for two data sets and models. The Virgo 
spectra (red symbols with error bars), the SDSS spectra (green symbols with error bars), 
single stellar population models (blue grid lines) from bottom to top of ages 17, 12, 8, 5, 3, 2, 
and 1.5 Gyrs and models from right to left of metallicities 0.5, 0.25, 0, —0.25, —0.5, —1, —1.5 
and —2 (pink grid lines) are plotted. The fits to the index values for the SDSS galaxies after 
the correction for hydrogen emission are shown as a light blue line in each panel. 
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The fits for H7 and B.S do not look as good wlien compared to tlie Virgo measurements 
(see Figure 3.1). Tlie reason for tliis is most likely the sensitivity of the H'y and H6 indices 
to changes in individual elemental abundances, with perhaps a non-negligible contribution 
from systematic errors between the two data sets. Elements N, C, and O have a profound 
and interdependent effect on this spectral region, and suggest themselves as candidates for 
further investigation. 



6. Discussion, Summary, and Conclusion 



Using the correction factor {-ft-) as derived above and associating it with the linear offset 



in Ha this work has shown that for SDSS spectra cleaned like those in lGraves et al.l ( 120071 ). 
a emission correction factor on the order of 0.5 A in Ha and 0.2 A in H/9 needs to be 
applied to the quiescent galaxies in order to better determine the mean age and metallicity 
of these galactic av erages. These correction factors are much larger than those estimated in 



Graves et al.l (120071 ). where the estimated correction factors from measurements of Oil give 
Ha = 0.082 A and H/3 = 0.027 A. This work also shows that the higher-order Balmer lines 
H7 and B.6 have such modest emission corrections that other uncertainties dominate the 
error budget. At least for our comparison sample of Virgo cluster elliptical galaxies, H7 
and B.6 suffer from contamination from the varying of other elemental abundances, making 
age determination from these indices more complicated. This makes the need for correct 
Ha and H/3 measurements all the more important, since they are relatively insensitive to 
changes in metal abundances. 



These results are in qualitative agreement with the results from lEisenstein et al.l (120031 ) 
who showed that in their SDSS spectra H/3 suffers from interstellar emission lines and 
speculated that non-solar abundance ratios were to blame for the differences in age 
determinations made from H/3, H7 and US. This is in agreement with these findings, where 
not only is there obvious hydrogen emission, but also it is clear that relative abundance 
ratios would certainly have to be taken into account in order to determine galactic ages 
from H7 and B.S. 

For this work the Ha correction was chosen as the basis for the corrections due to the facts 
that Ha is ~ 3 times more sensitive to hydrogen emission than H/3 and that the line fit for 
Ha as measured in the Virgo galaxies had a much tighter fit (RMS = 0.066 A) compared 
with H/3 (RMS = 0.15 A) . This lead to a corrected SDSS line fit with much less 
uncertainty, especially for galaxies inside the range 3.0 < Mg b < 4.3 A where the 
corrections for H/3 yield the most reliable results. Another reason for choosing Ha instead 
of directly using the correction that one could get from the H/3 plot is that we wanted to 
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preserve for age determination. 



The plot for H/3 in figure 3.1 still shows a discrepancy between the corrected SDSS data 
(light blue line) and the Virgo data (red line). Although, statistically speaking, the 
discrepancy is of marginal significance, this residual difference could be due to a few 
variables such as varying abundance ratios within the Virgo galaxies, slightly different 
decrements values than the ones used here, or the possibility that the Ha emission 
correction may be in part due to an difference in the mean ages of the two samples. 
Thomas et al.l (120051 ) showed that cluster galaxies tend to be around 2 Gyr older than field 
galaxies. 

A contributing factor might be that Ha and H/3 have different age sensitivities. In order to 
investigate the age sensiti vities, the Z ver s us ag e sensitivity parameter (Zsp) was calculated 
for both Ha and H/3 as in lWorthey et al.l (119941 ). The Zsp parameter is the modeled change 
in index strength due to a change in fractional metallicity (Z = 0.01689 x 10 

[Fe/H]^ 

divided by the change in index strength due to a change in fractional population age. 



Zsp 



[^U/(^log(Z))] 
[5L/((51og(age)] 



(6) 



Here 6lm/ Slog{Z) is the partial derivative of the index with respect to metallicity at age = 
12 Gyrs. Similarly, 5Ia/51og(age) is the partial derivative of the index with respect to age 
at solar metallicity. T hese sensitivities are shown in Table 3.3 along with the original 
Worthey et al.l ( 1l994j ) H/9 sensitivity. Note that the models indicate that both Ha and H/3 
are age indicators of the same sensitivity. This would mean that the difference between Ha 
and H/3 is more likely to be due to abundance ratios or the decrement values used than it is 
the age sensitivities of Ha and H/3. H owever, since both Ha and H/3 are insensitive to 
changes in abundance ratio variation (IServen et al.ll2005l ) and decrement values can change 
from the effects of local dust, the difference would most likely stem from the choice of 
decrement. However, the measured extragalactic decrement values tend to be larger than 
those used here and that using these larger decrements would only increase the observed 
residual difference! We are therefore forced to postulate that at least one other systematic 
offset may be operating between the Virgo and SDSS data sets. 
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Table 3.3 



Index 


Zsp 


Ha 


0.8 


H/3 


0.8 


H/3 (Worthey et al. 1994) 


0.6 



The applicability of this work for other grand SDSS averages may be limited due to the 
details of the sample selection. It is also unlikely to be applicable to individual red 
sequence galaxies due to the wide dispersion in index values and p ossible age ef f ects. One 



possible avenue could be to scale other averages to match those of iGraves et al.l (120071 ) . but 
it would be a better idea to echo the work shown here by comparing those averages with 
the minimal-emission Virgo data set and determining a new correction. 

It is conceptually possible to solve for both age and emission correction by considering both 
H/9 and Ha simultaneously, and increasing the emission correction until both indices give 
similar ages against a model grid. The obvious trouble with that scheme is that the solution 
then becomes model dependent. There is also a strong anticorrelation between derived age 
and emission correction. There is also the contribution of N emission near Ha, whose 
presence may cause spuriously large Ha emission measurements in individual galaxies. 

Speculation aside, it is safe to say that there is emission contamination in the SDSS spectra 
and that it is reasonably well accounted for by the linear fits presented in this work. In 
data sets that include Ha, observed OH and OIII do not need to be used as a proxy for 
Balmer emission. Emission contamination is much less of a problem for H7 and H5, but 
interpretation of these indices is complicated by the probable effects of individual elemental 
abundances. 

We would like to thank Genevieve J. Graves for providing the Sloan Digital Sky Survey 
spectra, as well as her advice and input for this paper. Major funding for this work was 
provided by National Science Foundation grants 0307487 and 0346347. 
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CHAPTER FOUR 
NH and Mg Index Trends in Elliptical Galaxies 



Introduction 



One of the biggest problems with any attempt to determine the chemical make up of a 
stellar system is trying to disentangle the effects of C, N, and O in any given spectrum. 
This is due to the fact that in cool stars and elliptical galaxies, C, N, and O show 
themselves in any spectrum almost exclusively through molecular species such as NH, CN, 
C2, CH, and CO rather than atomic species. The intertwining of these three elements 
starts with CO. The fact that CO has the highest dissociation energy of these molecules 
means that CO will form the most prolifically if given enough O and C. Since O is typically 
the most abundant of these three elements, C becomes incorporated into CO instead of the 
other molecular species. The rest of these molecules are connected through balancing of 
molecular equilibria. These interactions give a net effect of O acting like anti-CN since 
adding any O will d ecrease the amount of C available for the formation of other molecules 
JServen et ahlboosh . 



To begin disentangling these three elements it should be possible to start with pseudo 
equivalent width indic es that are sensitive to these elements such as C24668 



C04685 f Serven et al. 



(IWorthey et al.lll994b ) , whi ch is sensitive to C abundance, and the indices C05161 and 



2OO5I ) which are insensitive to N but react to C and O abundaces. 



After getting a grip on the C and O abundances, one can use the CN band to determine 
the N ab undance. Unfo rtunately, disentangling C, N, and O this way turns out to be 
difficult ( lBursteinll2003l ). To low precision, most previous results agree that Mg, C, N, and 
Na appear to be enhanced i n large elliptical galaxies and also correlated with velocity 



dispe rsion (IWortheylll998al : ISanchez-Blazquez et al.ll2003l : iKelson et al.ll2006l : iGraves et al. 

mm- 



What is most often suggested as an alternative for determining the N abundance is the use 
of the NH feature at 3360 A since, as has been noted before, it is insensitive to C and O 



(jSnederull973l: iNorris et al.l 120021) and it is also directly and se nsitively measuring N 
abundances (iBessell fc Norrid Il982l : iTomkin fc LambertI Il984j ) . 



Below, I will show that this is not the whole story and that there are other contributors to 
the NH feature. Unfortunately, until recently there had been fewer than 15 early-type 



This preprint was prepared with the AAS IAT5t;X macros v5.2. 
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galaxies with published N H3360 values due in large part to relative insensitivity of 
detectors in the near-UV (IToloba et al.ll2009l ). So making use of the NH3360 feature was 
almost impossible u ntil the introduction of NH3360 values of 35 early-type galaxies from 



Toloba et all (|2009|). 



In 



Toloba et al.l (120091) this sample of 35 galaxies was me asured using indice s NH3360 



f Pavidge fc Clarklll994h . CN03862, CN04175, C04685 flServen et al.ll2005h . and Mg b 
(iBurstein et al.lll984j ). Their findings included that there exists a flat relation between the 
NH3360 index and velocity dispersion. This seems to indicate that there does not exist a 
velocity dispersion relation for nitrogen, contrary to work done previously. For example, 
since the CN relation is stronger and tighter than the C24 668 relation, it seem s that there 
should be a positive [N/Fe] trend with velocity dispersion ( iTrager et al.lll998l ). 

It is toward attempting to explain this disparity that the rest of this chapter is aimed. In 
the first section, t he response o f the N H3360 feature to various elements is calculated using 
models as done in IServen et al.l (120051 ) but with better models. Then, after seeing Mg 
contamination in the index, two new indices are defined, one for the contaminate Mg, and 
the other for N. The responses for these new indices are then calculated. In §3 the NH3360 
index along with NH3375, Mg3334, CNi and Mg b are plotted and compared to models to 
determine if the effects of an old metal-poor stellar population can explain the observed 
index trends. Lastly, the results and conclusions are discussed in §4. We find that the 
presence of an old metal-poor population does account for the observed index trends. 



8. Analysis 



The method for deterr nining the response of N H3360 as wel l as th e two new indices was 
similar to that used in lServen et al.l (120051 ). In IServen et al.l (120051 ) simple models of a 
galaxy spectrum were constructed, using a G dwarf and a K giant. These models varied in 
that there was one base model of solar metallicity and then 23 variations, each one with a 
particular elemental abundance doubled. Then, the ratios of these spectra were taken to 
find any spectral influence due to any particular element. 

The difference here is that the numbers tabulated come from measurin g these indices i n 
full single ste l lar po pulation models. These models are a version of the IWorthey Jl994h fc 



Trager et al.l (119981 ) models that use a grid of synthetic spectra in the optical (ILee et al. 
2OO9I ) . Age, overall metallicity Z, and 23 individual elemental abundances can be varied 



independently in the models, and spectra for single-bu rst (simp l e) stel lar populations 
produced. In this chapter plotted results are based on IWortheyl (119941 ) isochrones, but all 



25 



results are confirmed with Padova (iBertelli et al 



1994f ) evolution. 



The grid of synthetic spectra is complete enough to predict nearly arbitrary composition. 
For the remainder of this paper whenever the term model is used it is with reference to 
these models, as they are the backbone of the methods used. 

The responses of NH3360 and the two new indices NH3375 and Mg3334 can be found in 
Table 4.3. What can be seen is that the NH3360 index, while indeed insensitive to C and 
O, is sensitive to Mg and to a lesser degree Fe and Ni. The ant i- correlation with Mg is due 
to a small Mg absorption feature located in the blue continuum passband. In order to 
remove this Mg dependence, a new index was defined; index NH3375, with an altered blue 
pseudocontinuum definition that avoids this feature. The new index NH3375 response, 
which is also found in Table 4.3, shows a smaller dependence on Mg with some small 
dependence on Ti and Ni. Unfortunately, the sensitivity of NH3375 to Ti and Ni comes 
from features of these two elements that overlap the NH feature itself. This limits the 
degree to which these sensitivities can be removed from any index definition. 

For the sake of complementarity, a new index was defined for the small Mg feature as well: 
index Mg3334. The response of this index is surprisingly clean, showing only small 
sensitivity to O (see Table 4.3). The definitions of these new indices, as well as NH3360, 
can be found in Table 14.11 The index sensitivities to N and Mg are modeled in Figure [H 
where the first panel shows NH3360, NH3375, and Mg3334 as functions of N abundance 
and the second panel shows these indices as functions of the Mg abundance. 




-0.3-0.2-0.1 0.1 0.2 0.3 0.4 0.5 -0.3-0.2-0.1 0.1 0.2 0.3 0.4 0.5 



[N/Fe] [Mg/Fe] 

Fig. 4.1.— In the first panel NH3360 (red), NH3375 (green), Mg3334 (blue) against [N/Fe]. 
In the second panel NH3360 (red), NH3375 (green), and Mg3334 (blue) against [Mg/Fe] for 
the second panel. This figure illustrates the sensitivity of these indices to changes in the N 
and Mg abundance at fixed [Z/H] = [Fe/H] = 0.0 and age = 8 Gyr with [Mg/Fe] = 0.0 for 
the first panel and [N/Fe] = 0.0 for the second. 
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Table 4.1 



Index 



Blue Passband 



Index Passband 



Red Passband 



Reference 



NH3360 3332-3350 
NH3375 3342-3352 
Mg3334 3310-3320 

Mg b 5142.625-5161.375 



3350-3400 
3350-3400 
3328-3340 
5160.125-5192.625 



3415-3435 
3415-3435 
3342-3355 
5191.375-5206.375 



Davidere fc Clark (1994) 
This work 
This work 
Trager et al. (1998) 



Table 4.1: Index passband definitions in A. 



A Z versus age sensiti vity para r neter (Zsp) was calculated for NH3360, NH3375, Mg3334, 



and Mg b as it was in IWortheyi ( 119941 ). The Zsp is the ratio of the percentage change in age 
to the percentage change in Z ( Z 0.01689 x 10 P^/"'^]) of the index measured as shown 
below. 

^ [^U/^log(Z)] 
[5L/51og(age)] 

Here 6lm/ Slog{Z) is the partial derivative of the index with respect to metallicity at age = 
12 Gyrs. Similarly, 5Ia/51og(age) is the partial derivative of the index with respect to age 
at solar metallicity. These sensitivities are shown in Table 14.21 Note that the models 
indicate that both NH3360 and NH3375 are far more sensitive to age than they are 
metallicity, especially NH3360 which has almost no metallicity sensitivity in the metal rich 
regime, while the Mg indices are more sensitive to metallicity. 

Table 4.2 



Index 


Zsp 


NH3360 


0.2 


NH3375 


0.6 


Mg3334 


1.1 


Mg b 


1.7 



Table 4.2: Table 4.2 shows the Z versus age sensitivity parameter (Zsp), which gauges how 
changes in metallicity and age effect various indices. A large Zsp indicates a larger depen- 
dence on the overall metallicity then on age with 1.0 indicating that age and metallicity 
effect the index equally. 

A plot illustrating the sensitivity of these indices to metallicity and age can be found in 
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Figure 14.21 Note that the NH3360 index, although age sensitive, goes nearly fiat as a 
function of metallicity. The other indices exhibit behavior of increasing and plateauing out 
over time. 




Fig. 4.2. — NH3360, NH3375, Mg3334, and Mg b against age for various metallicities. 
The NH3360 index (top left), the NH3375 index (bottom left), the Mg3334 index (top 
right), and the Mg b index (bottom right) are plotted for metallicities from top to bottom 
0.5, 0.25, 0, -0.25, -0.5, -1, -1.5, -2 and ages 1.5, 2, 3, 5, 8, 12 and 17 Gyrs. 



9. Observations and Results 



The iToloba et al.l ( 120091 ) sample consists of long-slit spectra for 35 elliptical galaxies 
collected with the 4.2 m William Herschel Telescope at Roque de los Muchachos 
Observatory, using the ISIS spectrograph. The spectra were extracted within a central 
equivalent aperture of 4" and broadened to a velocity dispersion of 300 km s~^. 
Measurements for galaxies with velocity dispersions greater than 300 km s~^ were corrected 
to a central velocity dispersion of 300 km s~^ using their own models. The wavelength 
coverage is from 3140 to 4040 A with a resolution of 2.3 A (FWHM) and a typical signal to 
noise oi S/N = 40 per A. These spectra were chosen as a subset of those presented in 



Sanchez-Blazquez et al.l (l2006al ) so that this near-UV data could be supplemented with 
optical data in the wavelength range 3500 to 5250 A. The galaxies in this set were also 
chosen to include field, Virgo and Coma cluster ellipticals, which cover a range of velocity 
dispersions (130 < a < 330 km s~^). 

From this sample the NH3360, NH3375, Mg3334, CNi and Mg b Lick style indices were 
measured. The measurements for NH3360, NH3375 and Mg3334 were taken from the 
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Toloba et al.l (120091) spectra , while the measurements for CNi and Mg b were taken form 



the ISanchez-Blazquez et aLl ( l2006al ) optical spectra. These measurements were then 
plotted against the Fe4383 index to compare the trends of these near-UV indices which are 
sensitive to N and Mg with the slightly redder CNi and Mg b which are also sensitive to N 
and Mg (Fig. 14. 3p . For comparison, Mg3334 vs Mg b is also plotted. Th e trend lines shown 



in Figure 4.3 are best-fit hues calculated using fitexy.f (jPress et al.lll992l ). a program for 
finding the best-fit line for data with errors in both the x and y coordinates. It minimizes 
the distance of each point form the line while taking into account weighting by the 
precision of the individual measurements in both the x and y coordinate. 

The index plotted in the first panel of Fig. 14.31 is NH3360, in the second NH3375 and in the 
third is Mg3334, all against the Fe4383 index. The index plotted in the fourth panel of Fig. 
14.31 is CNi, in the fifth is Mg b both against the Fe4383 index and in the sixth is a plot of 
Mg3334 vs Mg b. In each panel along with the index are plotted three lines. These lines 
represent the respective indices measured from a 12 Gyr galactic model of solar metallicity 
[Fe/H] = 0.0 to a metallicity of [Fe/H] = 0.25. The red line is a model of solar metallicity 
with no subpopulation. The green line is a model of solar metallicity, but with 5 percent of 
the galactic mass consisting of a 12 Gyr, metal-poor ([Fe/H] = —1.5) subpopulation. The 
blue line is the same model except that the subpopulation is now 10 percent of the total 
mass. The light blue points are galaxies suspected as having a spectral calibration issue, as 
they have a large CNi but a small Fe4383 or small CNi and large Fe4383, something not 
seen in normal elliptical galaxies. 

What can be seen in Fig. 14.31 is that, for all three near-UV indices the index trends look 
fiat in contrast to the panels 4 & 5 of Fig. 14.31 which show a definite trend with increasing 
metallicity. What can also be seen is that, for the galactic models as the percent mass of 
the subpopulation is increased (red to green to blue) the model index trends flatten out; 
twisting to come into fairly good agreement with the observed indices within a metal-poor 
subpopulation fraction of 10 percent by mass. The redder indices of CNi and Mg b show 
that the subpopulation has little to no effect on these index trends except lowering the 
average metallicity. 

Further evidence for the existence of a old metal-poor population can be see in panel 6 of 
Fig. 14.31 what can b e seen here is that Mg3334 and Mg b which are fairly clean indicators 



of Mg ( IServen et al.l (j2005f ) & Table 4.3) do not appear to be measuring the same 
abundance trend. This is an indication that the near-UV indices do suffer from dilution 
due to an underlying bright and weak-lined near-UV population. 
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Fig. 4.3. — In the first, second and third panels, NH3360, NH3375 and Mg3334 indices vs 
Fe4383 are plotted respectively. In the fourth and fifth panels, CNi and Mg b indices vs 
Fe4383 are plotted respectively. In the sixth panel Mg3334 vs Mg b is plotted. In each panel 
along with the index are plotted three lines. These lines represent the respective indices 
measured from a 12 Gyr galaxy model of metallicity [Fe/H] = 0.0 to a metallicity of [Fe/H] 
= 0.25. The red line is a model of solar metallicity with no subpopulation. The green line 
is a model of solar metallicity, but with 5 percent of the galaxy mass consisting of a 12 Gyr, 
metal-poor ([Fe/H] = —1.5) subpopulation. The blue hne is the same model except that 
the subpopulation is now 10 percent of the total mass. The light blue points are galaxies 
suspected as having a spectral calibration issue, as they have a large CNi but a small Fe4383 
or small CNi and large Fe4383; something not seen in normal elliptical galaxies. 



10. Discussion, Summary, and Conclusion 

With the use of simple stellar population models it has been demonstrated have that the 
NH3360 index, while being C and O insensitive, suffers from contamination from other 
elements, most notably Mg (see Table 4.3). This dependence of the NH3360 index on Mg 
can be effectively removed by a simple redefinition of the NH3360 index blue continuum 
passband to form a new index, index NH3375. This new index is now not only insensitive 
to C and O but also Mg (Table 4.3) and even though contaminations from other elements 
exist (e.g.,Ti), NH3375 is a cleaner measure of N than NH3360. 

The predicted behavior of NH3375, as seen in FiglH is an increased sensitivity to N 
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abundance (Table 4.3). This prediction is borne out in the data (see Fig. 14. 3p as an 
increase in the magnitude of the NH3375 measurements compared too that of NH3360. 
Also borne out in the data is that the index trends of NH3375 and NH3360 are nearly flat 
as is the Mg3334 index. This flat behavior is in contrast to the index trends seen in CNi 
and Mg b, the premise being that if these indices are measuring N and Mg, respectively, 
then these indices should show similar behaviors. 

There appears to be a difference between the underlying stellar populations that produce 
the near-UV spectra and the visible spectra at least as regards these near-UV indices. This 
underlying near-UV spec trum is likely due to a fraction of metal-poor population on the 
order of 5-10 % by mass (IWorthey et al.lll996l ). The addition of a metal-poor 
subpopulation to our models does account for the discrepancies observed in the index 
trends (Fig. 14. 3p most notably that of M g3334 and Mg b . The low metal-poor fractions 
tested here are of order those inferred by IWorthey et al.l (119961 ) for the amount of mass 
locked into low metallicity stars in present day galaxies. 

The Mg3334 index turns out to be a useful tool as it is a fairly clean index with very little 
contamination (Table 4.3), which allows for a direct comparison to Mg b (Fig l4.3p . The 
remarkable fiat behavior of Mg3334 compared with Mg b can therefore not be attributed to 
Mg abundance, but therefore must be due to some other effect, and the metal-poor 
subpopulations hypothesis explains the difference nicely. Indeed, the pair of the indices 
taken together can be used to characterize the amount and abundance spread in the 
underlying metal poor fraction of stellar mass. 



With the inherently composite nature of the near-UV spectra (IBurstein et al.lll988l ). with 
metal-poor light mingling with metal-rich light, deriving a N abundance from NH3375 or 
NH3360 looks more difficult than heretofore suspected. It is outside the scope of this 
chapter, but by using of Mg3334 it may be possible to uncover N abundance by comparing 
derived Mg abundances from Mg3334 and Mg b and then calibrating the underlying 
metal-poor population until the Mg3334 abundance agrees with that of Mg b. This 
metal-poor calibration could then be applied to NH3375, which in turn could be used to 
derive a N abundance for comparison with N abundances derived form C, N and O indices. 
Hopefully these measurements will be in agr eement with the N enhancements in la rge 
elliptical galaxi es deduced by other authors (iGraves et al.l 120071 : iKelson et al.ll2006 : 
Worthevlll998ah . 



The existence of a small scatter Mg - a (velocity dispersion) relation among large elliptical 
galaxies implies increasingly effective Mg enrichment from Type II supernovae (i.e. massive 
star explosions) with increasing galaxy size. This might be because the initial mass 
function is stronger in large ellipticals (causing more Type II supernovae and light element 
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enrichment), or a Type la delaye d timescale is a cause, or because stellar binarism depends 
on environment (jWortlieylll998bl ). The similarity of the N - a and Na - a relations to that 
of the Mg - a relation imply that the same mechanism (that of Type II supernovae) is 
responsible for the observed trends and that contributions from nucleosyathesis on the 
asymptotic giant branch do not contribute significantly. That is to say, N appears to be 
"prim ary" rather than "secondary" according to this line of reasoning ( iHenry fc Worthey 



19991 ). The C - a relation, on the other hand, see ms to be inte r media te between Mg, N, Na 



and Fe-peak elements (from Type la supernovae; iTrager et al.l ( 119981 )). This implies that 
the C abunda nce may come frorn cont ributions from both supernovae and the asymptotic 
giant branch (IHenry fc Wortheylll999l ). 



Major funding for this work was provided by National Science Foundation grants 0307487 
and 0346347. 



REFERENCES 

Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F., & Nasi, E. 1994, A&AS, 106, 275 
Bessell, M. S., & Norris, J. 1982, ApJ, 263, L29 

Burstein, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 1984, ApJ, 287, 586 

Burstein, D., Bertola, F., Buson, L. M., Faber, S. M., & Lauer, T. R. 1988, ApJ, 328, 440 

Burstein, D. 2003, Star Formation Through Time, 297, 253 

Davidge, T. J., & Clark, C. C. 1994, AJ, 107, 946 

Gonzalez, J. J. 1993, Ph.D. Thesis, 

Graves, G. J., Faber, S. M., Schiavon, R. R, & Yan, R. 2007, ApJ, 671, 243 
Henry, R. B. C., & Worthey, G. 1999, PASP, 111, 919 

Kelson, D. D., Illingworth, G. D., Franx, M., & van Dokkum, P. G. 2006, ApJ, 653, 159 

Lee, H. c, Worthey, G., Dotter, A., Chaboyer, B., Jevremovic, D., Baron, E., Briley, M. 
M., Ferguson, J. W., Coelho, P., & Trager, S. C. 2009, ApJ, submitted 

Norris, J. E., Ryan, S. G., Beers, T. C., Aoki, W., & Ando, H. 2002, ApJ, 569, L107 



32 



Press, W. H., Teukolsky, S. A., Vetterling, W. T. and Flannery, B. P. Numerical Recipes in 
FORTRAN. Cambrige University Press, Cambridge, 2nd edition, 1992. 

Sanchez-Blazquez, P., Gorgas, J., Cardiel, N., Cenarro, J., & Gonzalez, J. J. 2003, ApJ, 
590, L91 

Sanchez-Blazquez, P., Gorgas, J., Cardiel, N., & Gonzalez, J. J. 2006a, A&A, 457, 787 

Sanchez-Blazquez, P., et al. 2006b, MNRAS, 371, 703 

Serven, J., Worthey, G., & Briley M. M. 2005, ApJ, 627, 754 

Sneden, C. 1973, ApJ, 184, 839 

Toloba, E., Sanchez-Blazquez, P., Gorgas, J., & Gibson, B. K. 2009, ApJ, 691, L95 
Tomkin, J., & Lambert, D. L. 1984, ApJ, 279, 220 

Trager, S. C, Worthey G., Faber, S. M., Burstein, D., & Gonzalez, J. J. 1998, ApJS, 116, 1 
Valdes, F., Gupta, R., Rose, J. A., Singh, H. P., & Bell, D. J. 2004, ApJS, 152, 251 
Worthey, G. 1994a, ApJS, 95, 107 

Worthey, G., Faber, S. M., Gonzalez, J. J., & Burstein, D. 1994b, ApJS, 94, 687 
Worthey, G., Dorman, B., & Jones, L. A. 1996, AJ, 112, 948 

Worthey, G. 1998a, PASP, 110, 888 

Worthey, G. 1998b, Abundance Profiles: Diagnostic Tools for Galaxy History, 147, 13 



This preprint was prepared with the A AS macros v5.2. 



33 



Table 4.3 





lo 


a 
























Index 


(A) 


(A) 


C 


N 


O 


Na 


Mg 


Al 


Si 


S 


K 


Ca 


Sc 


NH3360 


3.221 


0.285 


-0.05 


5.18 


-0.52 


-0.03 


-3.43 


-0.04 


-0.30 


-0.05 


0.00 


-0.20 


0.13 


NH3375 


4.688 


0.322 


-0.20 


4.51 


-0.87 


-0.03 


-0.71 


-0.04 


-0.04 


-0.05 


0.00 


-0.50 


-0.02 


Mg3334 


1.542 


0.105 


-0.43 


-0.22 


-1.14 


-0.04 


5.91 


-0.05 


0.14 


-0.04 


0.00 


-0.35 


-0.24 



Table 4.3: Table 4.3 shows the response of the NH3360, NH3375, and Mg3334 index defini- 
tions to various elements. The first column is the name of the index, the second column gives 
the index measurements in angstroms of equivalent width, and the third column gives the 
error associated with S/N = 100 at 5000 A. The remainder of the columns are the changes 
(enhanced minus unenhanced) in the index brought on by an element enhancement of 0.3 
dex ( or 0.15 dex for C) in units of the error of the third column. 



Table 4.3 Continued 



Index 


Ti 


V 


Or 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


Sr 


Ba 


Eu 


upX2 


NH3360 

NH3375 
Mg3334 


0.04 

-2.09 
-0.26 


0.37 

0.25 
-0.52 


-0.18 

-0.08 
0.18 


-0.12 

-0.25 
-0.65 


-1.70 

-0.78 
-0.29 


-0.41 

-0.27 
0.03 


1.36 

1.82 
-0.92 


-0.04 
0.00 
0.06 


-0.37 

-0.57 
-0.26 


0.01 

0.00 
-0.01 


0.00 

0.00 
0.00 


0.00 

0.01 
0.00 


-3.82 
-3.97 
3.10 
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CHAPTER FIVE 
The Effects of Velocity Dispersion on Index Trends 



11. Introduction 



As stated in the Chapter 1, IWorthey et al.l (1l992l ) presented models and observations 
indicating that individual elemental abundances for individual galaxies could be found by 
taking advantage of the fact that measured indices showed considerable scatter when 
compared to models. This scatter indicated that there were nonsolar abundance ratios 
within individual galaxies. 



Wortheyl (Il994j ) pointed out that this scatter along with the "3/2 rule" (Chap. 1 Sec. 1) 
for age and metallicity could be used to determine abundance ratio variations without 
having to know the age or metallicity to high precision if features could be found in the 
spectrum that were sensitive to such non-lockstep behavior. 



In IServen et al.l ( 120051 ) the task of finding non-lockstep features was undertaken with the 
use of synthetic spectra. The question of being to able to measure the effects of abundance 
ratio changes of 23 elements was explored with the use of element targeted indices, of those 
23, 18 looked measurable in spectra of integrated light spectra of S/N > 100. 

With this in mind, high quality spectra for 18 Virgo Clust er galaxies along wit h a Sloan 
Digital Sky Survey sample (Chapter 2) and a sarnple fromlToloba et al.l (120091) (Ch apter 3) 



were measured using the 25 l i ck ind ices (IWortheyl ll994J : IWorthey &: Ottaviani 



1997h. 1 R a 



index defined in I Cohen et al.l (119981 ) and 48 of the indices defined in IServen et al.l (120051 ) . 
These measurements were then analyzed to look for abundance trends with velocity and 
individual scatter in the abundance measurements as an indication of individual abundance 
variation. 

The rest of this chapter is organized as follows. The observations are discussed in §2, then 
the analysis of the spectra are outlined in §3. In §4, the results are presented and discussed. 



12. Observations 

The sample used for the measurements in this paper consist of new spectroscopic data for 
18 Virgo cluster galaxies, a set of averaged Sloan Digital Sky Survey spectra which have 
been coadded and binned by internal velocity dispersion into six stacks (Chapter 2) and a 
sample of the 35 Toloba spectra (Chapter 3). Both the Virgo and Sloan data were chosen 
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to be quiet, red galaxies, and without concern for orientation of the galaxies axis with 
respect to the instrument. 

The Virgo cluster galaxy spectra are high quality long slit spectra obtained with the Kitt 
Peak National Observatory (KPNO) Mayall 4m telescope with the Cassegrain 
spectrograph. These spectra were chosen primarily to cover a velocity dispersion from 50 to 
350 km/s. For each galaxy two overlapping spectrograph setups were used, one for the blue 
spectral range of 3200A to 5300A and one for the red spectral range of 5000A to 7500A. 

For the blue spectra the KPC-007 grating/grism with a 2 arcsec slit was used along with 
the T2KB CCD. For each galaxy, multiple exposures were taken. In the blue the exposure 
time calculator gave S/N = 40 for a 600 second exposure time for a galaxy at the faint end 
of the velocity dispersion range, so 6 to 8 such exposures were taken in the blue in an 
attempt to reach S/N = 100 after stacking the individual exposures. The time and number 
of exposures varied from galaxy to galaxy, but not by much. 

For the red spectra the Bl-420 grating/grism with a 2 arcsec slit was used along with the 
T2KB CCD. For each galaxy multiple exposures were taken. In the red the setup is more 
sensitive, and also the galaxies emit more red light so the 6 to 8 exposures were usually 
only 300 seconds in length, again with the idea to reach S/N =100 in the end. 

These spectra where then reduced using standard long slit spectroscopy reduction 
techniques with the use of IRAfB This produced one dimensional spectra binned linearly 
in wavelength. 

The SDSS spectra are high signal to noise spectra and are made up quiescent galaxy 
stacked spectra in 6 velocity dispersion bins. The bins are 70 < o" < 120 km s~^, 120 < o" < 
145 km s-\ 145 < a < 165 km s-\ 165 < a < 190 km s-\ 190 < a < 220 km s~^ and 220 
< a < 300 km s~^. In each bin are several hundred stacked spectra of similar galaxies. To 
read more about se lection criteria and image processing for the SDSS spectra see 
Graves et al.l koOl\ ]. 



The iToloba et al.l ( 120091 ) spectra were chosen to include field, Virgo and Coma ellipticals, 
spanning a range of velocity dispersions (130 < a < 330 km s~^). The index measurements 
come from two sources. For indices bluer than 3900 A the measurements were taken from 



spectra presented in iToloba et al.l (120091). The rernaining Toloba index measurements come 
from spectra presented in ISanchez-Blazquez et al.l ( l2006l ). The same set of galaxies is 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under Cooperative agreement with the National 
Science Foundation. 
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represented in all index measurements. 



13. Analysis 



For the purpose of comparing trends in the SDSS, Toloba and Virgo spectra, the spectra 
where broadened to a velocity dispersio n of 300 km s ~^. Measurements of the Lick indices 



and the Lick style indices introduced in IServen et al.l ( 120051 ) where then taken 



These measurements where then plotted first against the internal velocity dispersion of the 
individual galaxies. For the Virgo spectra, velocity dispersions and errors from 
Davies et. al.l (119871 ) were used, while center values of the velocity dispersion bins and half 
of the bins interv al was used for the v elocity dispersion and errors of the SDSS data, 
respectively. The iToloba et al.l ( l2009l ) values were calculated f rom the spectra using 



templ ates and the MOVEL and OPTEMA fitting algorithms (ISanchez-Blazquez et al. 
(120061) and references therein). 



Next these measurements and their errors were plotted against the Mg b and <Fe> indices. 
The first set should reveal correlations with alpha element like behavior and also with 
velocity dispersion, given the positive correlation between Mg b and velocity dispersion. 
The second set should reveal correlations of these indices with Fe. Note, however, that 
young age will push galaxies to weak metal line strength, regardless of element. 

It is also important to note that for errors in the large index measurements of the Virgo 
sample a error floor was imposed. This error floor was imposed due to the fact that for 
large indices the errors tend to average out giving unrealistically small errors. The error 
floor was constructed by fitting a line to the index measurements vs. index size, and then 
correcting those indices that fell below this line. This was done for CNi, CN2, Mgi, Mg2, 
TiOi, TIO2 (corrections in magnitudes ~ 0.03) and CaHK (corrections in angstroms ~ 0.3 

A). 

In these plots a tr end line was then fit in order to characterize the correlations. This was 
done with fitexy.f (jPress et al.lll992l ) which is a fortran program for finding the best fit line 
from data with errors in the x and y coordinate. Where the data for the two data sets is in 
agreement a single trend line is used, but there are indices for which two trend lines must 
be used due to differences in the data sets. Differences of note are the Balmer series, which 
is most likely due to hydrogen emission, TiO and CN data shifts most likely due to 
aperture differences or spectral shape errors. Aperture differences come from the fact that 
the Virgo spectra were collected using a narrow slit and thus contain light from the galactic 
nucleus. The SDSS spectra on the other hand contain more light from non-nuclear regions 
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of the galaxies. The spectral shape errors may come from differences in reduction 
techniques. While a distinct worry, this only seems to effect a few of the 74 indices. 

Plotted alongside everything else are 5 of the Virgo cluster galaxies with definite hydrogen 
emission. These can be seen as the gray crosses in the Appendix B figures. Also i n the 
graph s of index vs Mg b or <Fe> a model grid is plotted. These models are from IWorthey 
(]l994j ) plus my own additions, and span a range of age from 1.5 to 17 Gyr and metallicity 
from [Fe/H] =— 2.00 to 0.50. The 74 index plots can be found in Appendix B. 



14. Discussion, Summary, and conclusion 

What can be seen in every plot of Appendix B is that there is always a scatter in the data. 
This is a clear indication that there are individual elemental abundances for e ach galax y 



and is to be expected as it indicates the non-lockstep behavior pointed out in IWorthey 



(ll9M) 



Also evident from the plots is that there is by no means an overall consistency in the a, Mg 
b, or Fe trends. Even indices designed to target the same element can have greatly different 
trends, such is the case for Mg (App. B Figures 16 - 20), Ca (App. B Figures 26 - 28) and 
Ni (App. B Figures 60 - 65). Since these index sets generally range from blue wavelengths 
below 4100 A to much redder wavelengths, a probable explanation is that of the old metal 
poor population discussed in Chapter 3, which would have the effect of flattening out the 
bluer indices. A mere flattening, however, is not enough to explain the complete reversal in 
trends seen for Mg3835 (App. B Fig. 16) and CaHK (App. B Fig. 26). 

It is important too remember that often these indices are not pure indicators of a given 
element so inferring to much about a given element's behavior in general should be 
avoided. However, there are indices that are very clean and do give a fair representation of 
an element's behavior, such as Mg b, CaHK, Fe5406 and Na D. 

In the case of Ha and H/3 (App. B Figs. 1 and 2) the trends indicate that the relative 
intensity of the absorption feature decreases as a function of a, Mg and Fe. This general 
indicates that larger galaxies (i.e. larger velocity dispersion, a) tend to be older as 
hydrogen absorption decreases as stellar population grows older. This conclusion is 
supported by the Fe5406 index (Fig. 54), which indicates that the overall Fe content 
increases with increasing galaxy size. 

From Mg b it can be inferred that Mg tends to be enhanced with respect to Fe in larger 
older stellar populations. The Mg enhancement can be seen in Figure 20 as the index trend 
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in the third panel pulls up and away from the model grid. 



The CaHK index is of particular interest du e to its appare n t unde rabundance in larger 



galaxies. This effect has been seen before in iThomas et al.l (j2003bl ). who speculate that the 



increase of Ca underabundance with galaxy mass balances the higher total metallicities of 
more massive galaxies, so that calcium abundance is constant and does not increase with 
increasing galaxy mass. They argue that metallicity dependent supernova yields are the 
most promising explanation for the under abundance. 

Although there is a wealth of information that these indices hint at, full modeling of these 
galaxies is required to decouple the various elemental influences present in each index before 
anything definite can be said about the individual elemental abundances. The modeling of 
these galaxies through the use of these indices will be the subject of future work. 
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CHAPTER SIX 
Line Strength Gradients in Virgo Cluster Galaxies 



15. Introduction 



It has been well esta blished that l ine-strength profiles in elliptical galaxies vary with 
radius. Starti ng with Faber IQTsh the number of line strength gradients has grow n 



considerably ( DavidgJ 



1992 



Davies et al. 1993; Fisher et al. 1995; Gonzalez 1993 



Kuntschnerlll998l : iKuntschner et al.ll2006l ). The results of these investigations and others 
show some general results. 

One of these results is that Mg and Fe line-strength gradients decline with increasing 
radius although they are generally shallow. This is mostly interpreted as a decline in the 
overall metallicity and suggest that abundance ratios sta y more or less constant within a 
given elliptical galaxy (lGonzalezlll993l : iKuntschnerl Il998[ l . 



Results al so indicate that the H/j index, an age indicator, is generally con s tant a s a function 
of radius (lDavidgdll992l : lGonzalezlll993l ) with large scatter. iDavies et al.l (Il993l ) concluded 
that this implied no gradients in the age of the stellar population within individual 
ellipticals. That conclusion was contradicted by iHenry fc Wortheyl (jl999l ). who argued that 
the flatness of H/3 together with the fact that the metallic absorption features increase in 
strength toward the centers of galaxies implied that gala xy centers are rn ildly younger. 
Differences in the strengths of the H/3 measurements led iGonzalea (119931 ) to suggest that 
significant differences in ages of elliptical galaxies exist, even within a single galaxy. 

However, these results have been limited by the number of element sensitive indices and 
the chemical flexibility of the available population models to accurately determine the age 
and metallicity gradients of these galaxies. 

In this chapter the 18 Virgo cluster galaxies and 74 indices presented in Chapter 4 have 
been measured and analyzed to look for index trends within these galaxies as a function of 
radius. 

The rest of this chapter is organized as follows. The basic galaxy sample data is presented 
in §2, then the analysis of the spectra are outlined in §3. In §4 the results are presented 
and discussed. 
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16. Observations & Galaxy Data 



The sample used for the measurements in this chapter consist of the spectroscopic data 
presented in Chapter 4 for the 18 Virgo cluster galaxies, with the exception of processing 
the data to look for radial trends. The spectra were collected with the Kitt Peak National 
Observatory (KPNO) Mayall 4m telescope with the Cassegrain spectrograph. 

To look for radial trends the spectra were reduced using standard long slit spectroscopy 
reduction techniques with the use of IRAfI. This produced two dimensional spectra 
binned linearly in the wavelength and spacial directions. 

The basic galaxy parameters used can be found in Table 6.1 and Table 6.2. The data has 
been gathered from various sources. 

Table 6.1 has the following columns. 

Column 1: Galaxy catalog number. 

Column 2: Adopted effective radius (Rg) in arc seconds (calculated from De Vaucouleurs 



effective diameter (Ag) (IFaber et al.lll989l )). 
Column 3: Adopted effective radius (Rg) error in arc seconds (calculated from error 
estimates in 



Faber et all (119891)). 



Column 4: Velocity dispersion (cr) in km s ^ fiFaber et al.lll989l ) 



Column 5: Velocity dispersion (cr) error in km s~^ ( iFaber et al.l 119891 ). 
Column 6: Indicates if the galaxy shows hydrogen emission in Ha (this work). 
Table 6.2 has the following columns. 
Column 1: Galaxy catalog number. 
Column 2: Position angle (PA) in degrees. 
Column 3: Minor to major axis ratio (b/a). 

Column 4: Position angle and minor to major axis ratio references 



•^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under Cooperative agreement with the National 
Science Foundation. 
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Tal)le 6.1 



NGC 


Re 


Re Error 


a 


a Error 


emission 




arc seconds 


arc seconds 


km s~^ 


km s~^ 


in Ha 


1400 


37.77 


<3.3 


250 


25.0 


emission 


1407 


71.96 


<3.7 


285 


28.5 




2768 


53.35 


<3.7 


198 


27.7 


emission 


3156 


49.79 


<3.7 


112 


15.7 


emission 


3610 


12.80 


>3.7 


159 


22.3 




4278 


32.89 


<3.3 


266 


37.2 


emission 


4308 


1.11 


>3.7 


88 


12.3 




4365 


57.16 


<3.3 


248 


24.8 




4406 


90.60 


<3.3 


250 


25.0 




4434 


18.50 


<3.3 


117 


11.7 




4458 


26.74 


<3.3 


106 


10.6 




4472 


104.02 


<3.3 


287 


28.7 




4473 


24.95 


<3.3 


178 


17.8 




4478 


14.03 


<3.7 


149 


14.9 




4486B 


3.07 


>3.7 


161 


22.5 




4486 


104.02 


<3.3 


361 


36.1 


emission 


4489 


32.15 


<3.7 


49 


4.9 




4564 


21.73 


<3.3 


153 


15.3 





Table 6.1: Virgo Cluster galaxy parameters. 
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Table 6.2 



NGC 


PA 

(degrees) 


b/a 


reference 




1400 


35.0 





900 


Skrutskie et al. f2006 


) 


1407 


60.0 





950 


Jarrett et al. f2003) 




2768 


92.5 





460 


Jarrett et al. f2003) 




3156 


49.5 





571 


Adelman-McCarthv et al. 


^2008) 


3610 


134.2 





776 


Adelman-McCarthv et al. 


^2008) 


4278 


21.8 





925 


Adelman-McCarthv et al. 


^2008) 


4308 


40.6 





926 


Adelman-McCarthv et al. 


^2008) 


4365 


45.0 





740 


Jarrett et al. f2003) 




4406 


125.0 





670 


Jarrett et al. f2003) 




4434 


44.9 





958 


Adelman-McCarthv et al. 


^2008) 


4458 


11.3 





900 


Adelman-McCarthv et al. 


^20081 


4472 


172.1 





399 


Adelman-McCarthv et al. 


^2008) 


4473 


95.0 





540 


Jarrett et al. f2003) 




4478 


160.5 





942 


Adelman-McCarthv et al. 


^2008) 


4486B 


101.3 





926 


Adelman-McCarthv et al. 


^2008) 


4486 


152.5 





860 


Jarrett et al. f2003) 




4489 


4.4 





917 


Adelman-McCarthv et al. f2008l 


4564 


50.0 





440 


Skrutskie et al. f2006 


) 



Table 6.2: Virgo Cluster galaxy parameters. 



17. Analysis 



To analyze the radial properties of these 18 Virgo galaxies, the two dimensional spectra 
from the observations were further reduced by extracting them through 61 single pixel wide 
apertures, so that the 31st aperture was centered on the peak of the spectral flux. This 
produced 61 one dimensional spectra that sampled different radii symmetrically about the 
center of the galaxies. 



The one dimensional spectr a were then measured using the 25 Lick in dices (IWortheylll994 : 
Worthey &: Ottavianilll997n. 1 Ha index defined in ICohen et al.l (Il998[ ) and 48 of the indices 
defined in IServen et al.l (|2005[ ). These measurements where then analyzed to check that 
they where indeed symmetric and then to look for index trends within individual galaxies. 
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To look for symmetry about the center of the galaxies (aperture 31) index measurements 
from symmetric pairs (e.g., apertures 30 and 32) where compared and found to agree 
confirming that the spectra where in fact symmetric about the center. The index 
measurements of symmetric pairs were then averaged together, except for the center 
aperture and its nearest neighbors (e.g., apertures 30, 31 and 32) which where averaged 
together as a set. The end result is a set of 14 data points for each index per galaxy. The 
14 data points come as a result of cutting off data points with large errors. 

The radial increments come from the Cassegrain spectrograph's angular resolution. The 
angular resolution is .69 arcsec / pixel in the spatial direction. Applying this to the 
aperture extracted spectra means that aperture 31 is at radius of arcsec, aperture 32 is at 
a radius of 1.38 arcsec, etc. These intervals are what is represented in the figures of 
Appendix C where they have been transformed into the log of units of the effective radius. 
The exception is in the first data point where the average of the distance between the 
center aperture and its nearest neighbors has been used. 

In plots of Appendix C a tr end line was then fit in order to characterize the correlations. 



This was done with fitexy.f (jPress et al.l 119921 ) which is a fortran program for finding the 



best fit line from data with errors in the x and y coordinate. 



18. Discussion, Summary, and conclusion 



What can be seen in the plots of Appendix C is that index trends tend to vary 
considerably from galaxy to galaxy and from index to index making any sort of general 
statement about them difficult to make without a long list of exceptions. However, some 
trends can be seen. For example, the Fe Profiles (App. C Pages 153- 162) do tend to be 
shallow implying that the metallicity sta ys relatively constant throughout the galaxies in 
agreement with what others have found (iGonzalezI Il993l : iKuntschnerl Il998l ) . The most 
notable exception is that of NGC 4486 on page 48, where the inner most data point takes 
an exceptional dive. This behavior is not confined to this index and its presence in many 
others suggests the presence of a decoupled galaxy core. 

Also in agreement with what others have found is that the Mg to Fe ratio is more or less 
constant with radius. This is implied by the slope of the radial measurements in Figure 6.1, 
where Mg b vs. <Fe> is plotted for different radii of four representative galaxies. The 
central measurements are in the top right with decreasing radii toward the bottom left, and 
they parallel constant age model sequences, as would be expected for global Mg 
enrichment. This observation, also implies that a mechanism for the chemical enrichment 
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in these galaxies, most likely Type I I superno v ae, is the same throughout the galaxy. This 
agrees with the conclusions found in Iwortheyl ( 1998 ). who speculated that the environment 
may play a role in three key ways. One, the upper initial mass function is stronger in large 
ellipticals (more Type II supernovae and thus more light elements). Two, the Type la 
timescale is delayed in large ellipticals (less Fe peak enrichment), or that stellar binarism is 
keyed to environment. 
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Fig. 6.1. — Mg b and <Fe> measurements of four representative galaxies (black, red, green 
and light blue) for various radii are plotted. Single stellar population models (blue grid lines) 
from bottom to top of ages 17, 12, 8, 5, 3, 2, and 1.5 Gyrs and models from right to left 
of metallicities 0.5, 0.25, 0, —0.25, —0.5, —1, —1.5 and —2 (pink grid lines) are also plotted. 
The central measurement for each galaxy is in the top right with decreasing radii toward the 
bottom left. The constant slope indicates that the [Mg/Fe] remains constant as a function 
of radius. 



Indices that are sensitive to C, N and O (App. C Pages 114 - 120) tend to decline with 
increasing radius suggesting that these elements are enhanced in the core of these galaxies. 

H/3 and Ha are relatively fiat showing little to no change with increasing radius. This 
coupled with the fact that most of the metallic indicators decline with increasing radius 
suggest that these galaxies are mildly younger toward the center as found by 



Henry fc Wortheyi (119991 ). This is illustrated in Figure 6.2 where the same 4 representative 
galaxies have been plotted as in Figure 6.1. In Figure 6.2, B.(3 vs. [MgFe] is plotted for 
different radii for the four representative galaxies. The central measurements are in the top 
right with decreasing radii toward the bottom left, and they cross constant age model 
sequences (horizontal blue lines), as would be expected for galaxy centers which are slightly 
younger than their outskirts. 
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Fig. 6.2. — H/9 and [MgFe] measurements of four representative galaxies (black, red, green 
and light blue) for various radii are plotted. Single stellar population models (blue grid lines) 
from bottom to top of ages 17, 12, 8, 5, 3, 2, and 1.5 Gyrs and models from right to left 
of metallicities 0.5, 0.25, 0, —0.25, —0.5, —1, —1.5 and —2 (pink grid lines) are also plotted. 
The central measurement for each galaxy is in the top right with decreasing radii toward the 
bottom left. The crossing of constant age lines indicates that the centers of these galaxies 
are slightly younger than the outskirts. 




The Mg indices (App. C Pages 122 - 126) are of interest as the index trend seen by other 
authors (jGonzalezlll993l : iKuntschnerl Il998l ) is only clear in the Mg3835, Mg4780 and Mg b 
indices. The other large Mg indices (Mg 1 and Mg 2) do show a decline with radius but 
only after an initial bump up from the center. The fact that this behavior is seen in almost 
all of the galaxies and not in Mg b suggests that there is a wavelength dependent focal 
variation in the spectra that may effect this region of the spectrum for indices that have 
wavelength spans of 400 A or more. 

Although there is a wealth of information that these indices hint at, full modeling of these 
galaxies is required to decouple the various elemental influences present in each index before 
anything definite can be said about the individual elemental abundances. The modeling of 
these galaxies through the use of these indices will be the subject of future work. 
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CHAPTER SEVEN 
The Results 

Throughout the course of these investigations, several physical implications about the 
nature of elliptical galaxies and there evolutionary history have been discussed and implied. 
These implications range from age and metallically trends across galactic populations to 
considerations of individual galaxy structure. 

In Chapter 2, the possibility of determining the existence of individual elemental 
abundances was investigated. It was found that individual elemental abundances for 
elliptical galaxies were in fact, observable. The presence of some of these individual 
abundance patterns was found in Chapter 5 and in Appendix B, where it is observable as 
scatter in any of the many figures found there. The implication is that, for each galaxy, an 
individual evolutionary history maybe needed to explain each galaxy's chemical profile. 

Also in Chapter 5 and Appendix B, what is apparent from the figures is that along with 
individual chemical profiles there are systematic chemical trends that cross galactic 
populations. These trends (i.e., Mg b and CaHK) are found to correlate with galaxy size 
and suggests that environment plays a significant role in the evolution of elliptical galaxies. 
The evidence suggest that these trends are an artifact of supernovae enrichment that 
depends via mechanisms that are ambiguous on environment. 

In Chapter 3, the important side topic of hydrogen emission in the elliptical spectra from 
the Sloan Digital Sky Survey was investigated. The presence of emission contamination 
was determined by direct comparison of index measurements between the survey spectra 
and spectra of 18 Virgo Cluster galaxies. It was found that the age determinations of the 
survey averages were 2 to 3 Gyr too young. This is an important constraint for the 
evolutionary picture due to the fact that the initial formation of galaxies must take place 
sooner on a universal time scale if galaxies are older, and puts tighter limits on late epoch 
star formation. 

In Chapter 4, the presence of an old, metal-poor stellar subpopulation in elliptical galaxies 
was investigated and found to account for observed index trends, especially in the 
terrestrial ultraviolet. The presence of this subpopulation is important for two reasons. 
One it certainly needs to be accounted for in order to accurately model galactic properties 
such as age. Two, this blue population needs to be taken into account for an accurate 
determination of the N abundance from the NH feature. 

In Chapter 6, individual Virgo Cluster galaxies were spatially dissected and with the use of 
index measurements several properties of elliptical galaxy structure was deduced. The 
physical implications are twofold. One is that the chemical make up of individual galaxies 
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tend to be uniform form the center to the edge, with a few exceptions such as decoupled 
cores. Two is that the cores of these eUiptical galaxies tend to be younger than the 
outskirts, meaning they had star formation later or for a more prolonged period. 

The physical implications of this aid in the understanding of the history and evolution of 
elliptical galaxies. 
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APPENDIX A 
INDEX RESPONSE TABLE 
These 6 tables give the sensitivity of each index to a given element. The 1st column of the 
table gives the indices names, the 2nd gives the index measurements in angstroms of 
equivalent width, and the 3rd giving the error associated with S/N = 500 at 5000 A. The 
next 23 columns split between tables are the changes (enhanced minus unenhanced in units 
of angstroms) in the index due to the elemental enhancement 0.3 dex (0.15 dex for C) of 
the element labeling the column. The 24th column is the change brought about in the 
index by increasing all the alpha elements by a factor of 2 (upX2). 



Table 1 





lo 


a 






















Index 


(A) 


(A) 


C 


N 


O 


Na 


Mg 


Al 


Si 


S 


K 


Ca Sc 



CNi 


0.036 


0. 


,005 


3.69 


6.52 


-7.73 


-0.23 


-0.67 


-0.27 


2.40 


-0.40 


-0.02 


-1.38 


0.04 


CN2 


0.097 


0, 


,006 


3.19 


5.86 


-7.58 


-0.25 


-0.81 


-0.30 


3.19 


-0.40 


-0.02 


-1.44 


-0.02 


Ca4227 


2.394 


0. 


,084 


-3.55 


-2.41 


-1.24 


-0.05 


1.61 


0.09 


0.19 


-0.21 


0.00 


6.57 


-0.43 


G4300 


6.984 


0. 


,154 


4.12 


-0.05 


-3.85 


-0.43 


-0.90 


-0.14 


-0.81 


-0.16 


-0.01 


1.00 


0.23 


Fe4383 


7.526 


0. 


,222 


0.23 


-0.38 


-3.09 


-0.10 


-1.18 


-0.35 


-0.80 


-0.18 


-0.02 


-1.11 


0.13 


Ca4455 


1.826 


0. 


,127 


-0.39 


-0.19 


-1.39 


-0.04 


-0.07 


-0.03 


-0.11 


-0.10 


0.00 


0.23 


0.00 


Fc4531 


4.450 


0, 


190 


-0.40 


-0.19 


-1.75 


0.06 


-0.55 


-0.06 


-0.81 


-0.11 


0.00 


-0.25 


-0.01 


C2 4668 


4.251 


0, 


,315 


4.88 


-0.45 


-5.08 


-0.23 


-0.70 


-0.02 


-1.81 


-0.24 


0.00 


0.02 


-0.11 


H/3 


1.153 


0, 


135 


-0.04 


0.13 


1.55 


0.24 


-1.88 


0.09 


-0.68 


0.10 


0.01 


0.03 


-0.07 


Fe5015 


5.375 


0. 


,298 


-0.82 


-0.17 


-1.44 


0.15 


-1.72 


0.07 


-0.35 


-0.15 


0.00 


0.13 


-0.04 


Mgi 


0.139 


0. 


,003 


7.73 


-0.91 


-10.15 


-0.67 


13.39 


-0.88 


0.24 


-0.48 


-0.06 


-1.18 


0.03 


Mg2 


0.315 


0. 


,004 


2.07 


-1.05 


-9.93 


-0.51 


19.32 


-0.78 


0.02 


-0.56 


-0.05 


-0.68 


-0.27 


Mg b 


4.882 


0. 


,149 


-1.10 


-0.33 


-2.03 


-0.35 


10.00 


-0.21 


-0.15 


-0.14 


-0.01 


0.02 


-0.01 


Fc5270 


4.149 


0, 


162 


-0.06 


-0.11 


-2.98 


-0.14 


-0.68 


-0.17 


-0.08 


-0.13 


-0.01 


0.80 


-0.14 


Fc5335 


3.700 


0, 


184 


-0.31 


-0.27 


-2.41 


-0.13 


-0.55 


-0.19 


-0.09 


-0.14 


-0.01 


-0.14 


0.08 


Fc5406 


2.475 


0, 


142 


-0.21 


-0.17 


-2.07 


-0.10 


-0.32 


-0.13 


-0.07 


-0.13 


-0.01 


-0.17 


0.03 


Fe5709 


1.176 


0. 


,119 


-0.22 


-0.19 


-1.10 


-0.45 


0.08 


-0.04 


-0.04 


-0.08 


0.00 


-0.07 


0.15 


Fe5782 


0.817 


0. 


,115 


-0.18 


-0.16 


-1.25 


-0.02 


0.18 


-0.03 


0.03 


-0.09 


0.00 


-0.06 


0.00 


Na D 


4.602 


0. 


,143 


0.22 


-0.36 


-5.04 


7.99 


0.49 


-0.16 


0.34 


-0.26 


-0.01 


-0.36 


0.01 


TiOi 


0.031 


0. 


,004 


-1.67 


0.17 


0.42 


0.14 


0.31 


0.14 


-0.17 


-0.14 


0.00 


-0.03 


0.00 


Ti02 


0.068 


0, 


,003 


-3.15 


-0.03 


0.15 


0.21 


0.52 


0.36 


-0.58 


-0.36 


0.00 


-0.55 


0.33 


HSa 


-5.490 


0. 


187 


-0.35 


-0.23 


5.79 


0.22 


0.60 


0.36 


1.00 


0.39 


-0.03 


1.34 


-0.26 


H7A 


-9.781 


0. 


,207 


-2.59 


0.37 


6.33 


0.37 


1.93 


0.39 


0.59 


0.37 


0.02 


-0.41 


-0.14 


RSf 


-0.600 


0. 


,120 


0.47 


0.23 


2.01 


0.05 


-0.07 


0.11 


2.32 


0.18 


0.00 


1.28 


-0.05 


H7f 


-2.258 


0. 


,120 


-2.92 


0.21 


4.35 


0.28 


1.17 


0.27 


0.29 


0.24 


0.01 


-0.41 


-0.38 


Ha 


0.660 


0. 


,137 


1.70 


0.47 


-0.22 


-0.09 


-0.71 


-0.08 


-0.06 


0.15 


0.00 


-0.23 


-0.04 
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Table 1 Continued 



Index 


1 1 


V 


or 


ivin 


re 


y^o 


iNl 


L^u 


Zjll 


or 


rsa 




upAz 




0.46 


-1.42 


-1.31 


-0.69 


-0.48 


-0.42 


0.46 


-0.02 


0.00 


-0.04 


0.06 


-0.02 


-7.81 


CN2 


0.58 


-0.65 


-0.96 


-0.54 


-0.61 


-0.58 


0.42 


0.00 


0.00 


0.07 


0.05 


0.00 


-7.32 




0.69 


0.62 


-0.24 


0.19 


1.60 


0.10 


0.11 


0.00 


0.00 


-0.30 


0.00 


0.00 


5.11 




2.34 


0.05 


-0.64 


-0.03 


-2.44 


0.04 


0.31 


-0.02 


0.00 


0.00 


0.00 


0.00 


-2.88 


XT',, /I 


0.40 


0.76 


0.17 


-0.57 


6.97 


-0.01 


0.15 


0.00 


0.00 


0.00 


0.02 


0.00 


-6.73 


ua44oo 


0.42 


0.89 


1.18 


0.77 


-1.53 


0.14 


0.49 


0.00 


0.00 


0.00 


0.00 


0.00 


-1.28 


Fe4531 


2.18 


0.04 


U.74 


-0.18 


0.48 


0.23 


0.20 


U.U2 


0.00 


0.01 


0.76 


0.02 


-1.71 


C24668 


1.29 


-0.42 


-0.33 


-0.52 


0.26 


-0.33 


0.20 


0.04 


0.04 


-0.05 


0.01 


-0.02 


-6.14 


XT (3 


0.33 


-0.06 


-0.45 


-0.10 


-0.62 


-0.13 


0.87 


0.00 


0.00 


-0.02 


0.00 


0.00 


0.40 


Fe5015 


2.59 


0.04 


-0.21 


0.08 


1.02 


-0.01 


0.71 


0.00 


0.00 


-0.03 


0.00 


0.00 


-0.62 


Mgi 


-1.24 


-0.06 


-0.64 


0.09 


-5.09 


-0.36 


1.12 


0.21 


0.00 


0.00 


-0.15 


0.00 


-5.15 


Mg2 


0.78 


0.05 


-0.73 


0.02 


-3.95 


-0.46 


0.44 


0.05 


0.00 


0.00 


-0.07 


0.00 


2.80 


Mg 


0.32 


-0.06 


-2.20 


-0.18 


-1.55 


-0.07 


-0.02 


-0.07 


0.00 


0.00 


0.00 


0.00 


6.25 




0.27 


-0.20 


0.05 


0.21 


3.30 


0.31 


0.15 


-0.01 


0.00 


0.02 


0.00 


0.00 


-3.54 


Fe5335 


0.19 


0.01 


0.79 


0.11 


3.27 


0.38 


0.03 


0.00 


0.00 


0.00 


0.00 


0.00 


-3.94 


Fe5406 


-0.07 


0.07 


0.61 


0.26 


3.02 


-0.05 


0.02 


0.00 


0.00 


0.00 


-0.01 


0.00 


-3.49 


Fe5709 


0.26 


-0.13 


0.41 


-0.03 


0.91 


0.00 


0.21 


0.13 


0.00 


0.00 


0.00 


0.00 


-0.83 


Fe5782 


-0.17 


0.11 


1.67 


0.04 


-0.28 


0.00 


-0.07 


0.19 


0.00 


0.00 


0.00 


0.00 


-1.46 


Na D 


-0.63 


0.02 


0.00 


0.00 


-0.37 


-0.01 


0.08 


0.00 


0.00 


0.00 


0.00 


0.00 


-8.41 


TiOi 


3.44 


-0.31 


-0.14 


-0.06 


-0.22 


-0.03 


-0.25 


0.00 


0.00 


0.00 


0.03 


0.00 


3.92 


Ti02 


5.36 


0.58 


-0.09 


-0.09 


0.15 


-0.12 


-0.09 


0.00 


0.00 


0.00 


-0.03 


0.00 


4.45 


HSa 


-1.08 


2.11 


0.10 


-1.39 


-9.96 


0.46 


-0.34 


0.00 


0.00 


-0.57 


-0.06 


-0.05 


8.69 


H7A 


-2.86 


-0.50 


0.85 


-0.23 


-4.99 


-0.02 


-0.41 


0.00 


0.00 


0.00 


0.00 


0.00 


7.25 


HSf 


-1.15 


1.16 


-0.12 


-0.88 


-4.96 


-0.40 


-0.17 


0.00 


0.00 


-0.61 


-0.09 


-0.10 


4.34 


H7f 


-0.64 


-0.08 


0.98 


-0.11 


-2.98 


-0.02 


-0.27 


0.00 


0.00 


0.00 


0.02 


0.00 


5.66 


Ha 


-1.67 


-0.06 


-0.03 


-0.01 


-0.65 


0.04 


-0.18 


0.00 


0.00 


0.00 


-0.02 


-0.01 


-2.47 
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Table 1 



Index 


lo 

(A) 


a 

(A) 


C 


N 


O 


Na 


Mg 


Al 


Si 


S 


K 


Ca 


Sc 


C04685 


3.699 


0. 


.334 


6.45 


-1.14 


-6.15 


0.08 


-0.28 


-0.01 


-2.34 


-0.23 


0.00 


-0.77 


0.15 


C05161 


0.352 


0, 


.086 


4.24 


-0.09 


-3.42 


-0.01 


1.91 


-0.13 


-0.51 


-0.03 


0.00 


-0.23 


0.01 


CN03862 


10.767 


0, 


,196 


9.39 


15.11 


-18.54 


-0.10 


3.36 


-0.13 


-2.91 


-0.87 


-0.01 


-6.44 


-0.02 


CN04175 


3.109 


0. 


.286 


3.92 


7.09 


-8.59 


-0.20 


-1.22 


-0.33 


2.47 


-0.42 


-0.01 


-0.94 


-0.14 


Na8190 


1.834 


0. 


.190 


0.11 


-0.48 


-1.77 


1.94 


-0.26 


-0.14 


0.27 


0.08 


-0.01 


-0.11 


0.02 


Mg3835 


6.768 


0. 


.160 


-3.40 


-4.30 


2.29 


0.08 


6.27 


0.18 


0.41 


0.08 


0.01 


-0.09 


0.11 


Mg4780 


1.854 


0. 


.163 


-0.66 


-0.16 


-0.76 


-0.47 


4.35 


-0.25 


0.49 


-0.09 


-0.02 


-0.31 


-0.26 


A13953 


-2.718 


0. 


.240 


-0.85 


-0.29 


-1.52 


0.00 


-0.22 


5.07 


-0.44 


-0.10 


0.00 


4.15 


-0.06 


Si4101 


0.273 


0. 


.096 


0.38 


0.11 


1.83 


0.02 


-0.20 


0.06 


2.60 


0.17 


0.00 


1.05 


0.02 


Si4513 


-1.400 


0, 


,166 


0.35 


0.36 


-0.04 


0.05 


-0.50 


-0.14 


2.87 


0.00 


-0.01 


-0.58 


-0.04 


S4693 


0.342 


0. 


,113 


2.14 


-0.19 


-2.17 


-0.10 


0.73 


0.00 


-0.35 


0.06 


0.00 


0.26 


-0.17 


K4042 


0.415 


0. 


,085 


-0.14 


-0.02 


-0.14 


-0.06 


-0.32 


-0.13 


0.13 


-0.01 


0.18 


-0.55 


-0.35 


Cahk 


18.122 


0. 


,365 


-1.67 


-2.79 


0.43 


-0.24 


-6.60 


2.19 


-2.29 


-0.20 


-0.12 


25.97 


-0.38 


Ca8542 


3.196 


0. 


,174 


-0.28 


-0.24 


-1.39 


-0.16 


-1.34 


-0.19 


-0.56 


-0.11 


-0.01 


5.69 


0.00 


Ca8662 


2.339 


0. 


,149 


-0.08 


-0.29 


-1.22 


-0.09 


-1.40 


-0.19 


-0.57 


0.02 


-0.01 


5.18 


0.00 


Sc3613 


4.723 


0. 


,101 


-3.20 


-3.34 


0.86 


-0.01 


-0.55 


0.02 


0.73 


-0.04 


0.00 


-0.07 


1.11 


Sc4312 


0.553 


0. 


,116 


2.03 


0.11 


-0.79 


-0.04 


-1.25 


-0.15 


-0.03 


0.01 


-0.01 


0.77 


0.33 


Sc6292 


0.487 


0. 


,199 


-1.51 


-0.80 


1.49 


0.04 


0.53 


0.16 


0.13 


0.00 


0.00 


0.06 


0.35 


Ti4296 


5.974 


0. 


,111 


4.49 


-0.19 


-5.64 


-0.23 


-0.88 


-0.16 


-1.02 


-0.29 


-0.01 


1.36 


0.32 


Ti4533 


1.544 


n 


,uoo 


-0.23 


-0.14 


-1.80 


-0.07 


0.34 


-0.03 


0.26 


-0.11 


0.00 


0.21 


-0.03 


Ti5000 


2.184 


0, 


,209 


-1.40 


-0.23 


-0.87 


0.31 


-0.44 


0.06 


-0.21 


-0.18 


0.00 


-0.06 


-0.11 


V4112 


-6.730 


0, 


,287 


-1.11 


-0.89 


5.05 


0.18 


0.78 


0.29 


1.03 


0.30 


0.01 


1.00 


-0.18 


V4928 


0.535 


0, 


,101 


0.58 


0.02 


-0.71 


-0.04 


-0.34 


-0.04 


-0.08 


0.01 


0.00 


-0.05 


0.03 


V6604 


0.199 


0, 


,117 


-0.40 


-0.55 


0.36 


-0.03 


-0.02 


-0.06 


-0.04 


-0.01 


0.00 


-0.04 


0.16 


Cr3594 


1.577 


0, 


,124 


-1.77 


-1.47 


0.53 


0.05 


0.64 


0.07 


0.49 


-0.04 


0.00 


0.08 


0.07 


Cr4264 


3.251 


0. 


,156 


0.00 


-3.21 


-1.22 


-0.13 


0.44 


-0.13 


0.47 


-0.11 


-0.01 


-0.50 


-0.07 


Cr5206 


1.688 


0, 


,085 


-2.02 


-0.38 


-0.59 


-0.15 


2.03 


-0.15 


0.60 


-0.13 


-0.01 


-0.22 


-0.09 


Mn3794 


-6.864 


0, 


,233 


0.08 


-0.13 


3.40 


0.06 


1.39 


0.07 


-0.19 


0.24 


0.01 


-0.19 


0.00 


Mn4018 


-14.657 


0, 


,305 


0.18 


-1.83 


3.63 


0.13 


2.06 


-3.82 


-0.46 


0.27 


0.05 


-14.36 


0.59 


Mn4461 


0.494 


0, 


,152 


-3.64 


-0.29 


0.62 


0.07 


1.74 


0.32 


-1.64 


-0.04 


0.02 


1.61 


-0.14 


Mn4757 


0.345 


0, 


,141 


-2.41 


0.08 


1.16 


0.23 


1.21 


-0.06 


0.47 


0.06 


0.00 


-0.05 


0.02 


Fe4058 


6.046 


0, 


,165 


-1.33 


-1.25 


-1.96 


-0.17 


0.03 


-0.29 


-1.26 


-0.15 


0.09 


-1.67 


-0.15 


Fe4930 


3.160 


0, 


,168 


-0.21 


-0.20 


-1.94 


-0.06 


-0.44 


-0.08 


-0.33 


-0.12 


-0.01 


-0.11 


0.08 
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Table 1 Continued 



Index 


Ti 


V 


Cr 


Mn 


re 


Co 


IN 1 


Cu 




or 


ria 


IjU 


upAz 


CU4b85 


0.15 


-0.74 


0.59 


-0.32 


0.30 


0.01 


0.12 


0.03 


-0.01 


-0.03 


-0.13 


0.02 


-9.65 


C05161 


-1.85 


0.01 


-1.21 


-0.15 


-1.36 


0.06 


0.03 


0.05 


0.00 


0.01 


0.01 


0.00 


-4.45 


AT O O / ■ O 


-0.83 


0.59 


0.57 


-0.11 


-2.62 


0.40 


-0.93 


0.00 


0.00 


0.00 


0.00 


-0.08 


-25.83 


/""AT /I 1 '7r 

ClNU4i75 


-0.22 


-1.40 


-1.75 


-0.53 


0.88 


-0.59 


0.43 


-0.01 


0.00 


0.17 


0.08 


0.09 


-9.67 


Na819U 


-2.07 


0.02 


0.02 


0.01 


0.10 


1.37 


0.00 


-0.03 


0.00 


0.00 


0.00 


0.00 


-4.75 




-1.07 


0.10 


-0.73 


0.56 


0.06 


0.46 


-0.85 


0.00 


0.00 


0.00 


0.00 


0.01 


6.82 


Mg478U 


0.84 


-0.52 


0.23 


0.03 


0.28 


0.16 


-0.32 


0.00 


0.00 


-0.10 


0.00 


0.00 


2.73 




-0.24 


-0.24 


-0.35 


-0.49 


-2.78 


0.01 


-0.02 


0.00 


0.00 


0.00 


0.00 


-0.04 


0.15 


Si4101 


-1.22 


0.39 


-0.29 


-0.37 


-4.20 


-0.45 


-0.22 


0.00 


0.00 


-0.17 


-0.06 


-0.08 


3.98 


Si4513 


0.03 


-0.81 


-0.17 


-0.70 


0.64 


0.10 


-0.28 


0.02 


0.00 


0.01 


-0.50 


0.02 


1.13 


S4693 


0.11 


-0.23 


0.37 


-0.06 


-1.49 


0.03 


-0.08 


-0.02 


-0.05 


-0.02 


0.02 


-0.01 


-1.55 


K4U4z 


-1.09 


0.29 


0.12 


1.71 


2.80 


-0.58 


-0.56 


-0.13 


0.00 


-0.01 


0.00 


0.00 


-2.40 


Cahk 


2.25 


-0.30 


-0.12 


-7.49 


-9.56 


0.86 


-0.39 


-0.12 


0.00 


0.00 


0.03 


0.05 


18.17 


Ca8o42 


0.08 


0.00 


0.03 


0.00 


-1.57 


-0.15 


0.00 


0.00 


0.00 


0.00 


0.01 


0.00 


2.56 


CasDoJ 


-0.55 


0.00 


-0.03 


-0.03 


-1.03 


-0.01 


-0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


1.57 


Sc3613 


-0.83 


-0.75 


0.72 


-0.10 


6.90 


-0.45 


1.53 


0.00 


0.00 


0.00 


0.00 


0.00 


-0.22 


Sc4312 


-0.75 


-0.27 


-2.14 


0.09 


1.82 


-0.03 


-0.44 


0.00 


0.00 


0.02 


-0.01 


0.00 


-1.83 


bc6292 


0.81 


0.13 


-0.16 


-0.01 


0.21 


0.06 


-0.20 


0.00 


-0.03 


0.00 


-0.03 


0.01 


2.94 


i i4z9d 


3.15 


0.05 


-0.88 


-0.18 


-1.93 


0.03 


0.53 


-0.01 


0.00 


0.01 


0.00 


0.00 


-4.00 


Ti4533 


1.95 


0.14 


0.68 


-0.23 


1.01 


0.43 


-0.50 


0.04 


0.00 


0.01 


-1.05 


0.01 


-0.12 
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u.uu 
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l.Do 


V4112 


-0.95 


1.91 


0.42 


-0.73 


-8.03 


0.17 


-0.27 


0.00 


0.00 


-0.25 


0.07 


0.09 


7.62 


V4928 


-0.55 


0.13 


0.16 


-0.04 


-0.18 


-0.06 


0.10 


0.00 


0.00 


-0.04 


0.39 


0.00 


-1.67 


V6604 


-0.23 


0.13 


-0.10 


0.00 


0.79 


0.03 


0.18 


0.00 


0.00 


0.00 


0.00 


0.00 


0.26 


Cr3594 


0.31 


0.01 


4.49 


-0.15 


-3.69 


-0.41 


1.83 


0.00 


0.00 


0.00 


0.00 


-0.01 


1.42 


Cr4264 


0.35 


-0.26 


4.07 


0.48 


1.61 


0.04 


-0.30 


0.02 


0.00 


-0.01 


0.00 


0.00 


-1.75 


Cr5206 


1.28 


-0.05 


5.26 


-0.07 


-1.19 


-0.26 


-0.30 


0.03 


0.00 


-0.01 


0.00 


0.00 


0.81 


Mn3794 


-2.20 


0.10 


-0.07 


0.47 


-5.41 


-0.44 


0.93 


0.00 


0.00 


0.00 


0.00 


0.00 


3.16 


Mn4018 


1.67 


-1.19 


0.40 


4.72 


6.00 


0.16 


0.45 


0.11 


0.00 


0.00 


0.03 


0.00 


-5.06 


Mn4461 


0.82 


0.80 


0.50 


2.11 


-1.22 


-0.11 


-0.16 


-0.01 


0.00 


0.00 


-0.13 


-0.02 


3.13 


Mn4757 


-0.43 


0.09 


0.44 


1.36 


0.02 


0.21 


-0.12 


0.00 


-0.25 


0.00 


-0.02 


0.00 


1.80 


Fe4058 


-0.86 


-0.75 


-0.38 


0.95 


11.72 


-0.40 


-0.61 


-0.15 


0.00 


0.77 


0.00 


-0.02 


-6.37 


Fe4930 


0.71 


-0.01 


-0.01 


0.00 


2.36 


-0.14 


0.17 


0.00 


0.00 


0.00 


0.23 


0.00 


-2.36 
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Table 1 





7o 


(T 
























Index 


(A) 


(A) 


C 


N 


O 


Na 


Mg 


Al 


Si 


S 


K 


Ca 


Sc 


Co3701 


-0.231 


0.126 


-0.11 


-0.06 


-0.28 


0.00 


-0.08 


0.00 


-0.24 


-0.01 


0.00 


0.29 


-0.12 


Co3840 


3.279 


0.143 


-1.86 


-2.69 


1.68 


0.04 


4.55 


0.09 


-0.06 


0.08 


0.00 


-0.05 


0.11 


Co3876 


-0.530 


0.205 


3.28 


2.48 


-2.23 


-0.01 


0.79 


0.00 


-0.29 


-0.02 


0.00 


-3.78 


-0.03 


Co7815 


2.034 


0.272 


-0.65 


-0.49 


-0.37 


-0.06 


-0.44 


0.14 


-0.15 


-0.06 


-0.01 


-0.15 


0.00 


Co8185 


1.753 


0.199 


0.15 


-0.36 


-1.82 


1.83 


-0.31 


-0.12 


0.19 


0.05 


-0.01 


-0.07 


0.01 


Ni3667 


-5.338 


0.205 


1.09 


0.16 


2.50 


0.07 


1.04 


0.12 


-0.39 


0.15 


0.01 


-0.07 


0.11 


Ni3780 


-7.051 


0.255 


2.13 


1.03 


3.09 


0.01 


0.88 


0.04 


-0.42 


0.25 


0.00 


-0.16 


0.00 


Ni4292 


4.160 


0.132 


2.35 


-0.16 


-3.22 


-0.24 


-0.30 


-0.05 


-0.68 


-0.18 


0.00 


0.72 


-0.01 


Ni4910 


-0.650 


0.134 


0.23 


0.08 


0.60 


0.01 


-0.01 


0.01 


-0.14 


0.06 


0.00 


-0.07 


0.10 


Ni4976 


0.023 


0.123 


-0.06 


0.04 


0.19 


0.46 


1.02 


0.07 


0.26 


-0.03 


0.00 


0.02 


-0.08 


Ni5592 


0.996 


0.131 


0.27 


-0.27 


-1.14 


-0.07 


-0.58 


-0.07 


0.22 


0.01 


0.00 


2.51 


0.01 


Cu5217 


-2.304 


0.108 


-1.44 


0.04 


2.76 


0.20 


-6.18 


0.28 


0.12 


0.07 


0.01 


-0.42 


-0.02 


Cu5780 


0.674 


0.082 


-0.12 


0.00 


-1.52 


-0.02 


0.08 


-0.03 


0.08 


-0.11 


0.01 


-0.07 


-0.04 


Zn4720 


-0.635 


0.118 


1.97 


0.02 


-0.47 


-0.07 


-1.13 


0.06 


-0.39 


0.01 


0.00 


0.08 


0.20 


Ba4552 


0.023 


0.082 


-0.04 


0.11 


0.18 


0.08 


-0.65 


0.01 


-0.86 


0.00 


0.00 


-0.36 


0.01 


Ba4933 


-0.197 


0.117 


0.77 


0.09 


-0.16 


-0.15 


-0.54 


-0.02 


-0.03 


0.06 


0.00 


-0.01 


0.08 


Ba6142 


0.143 


0.117 


0.72 


0.01 


-0.22 


0.03 


-0.40 


-0.08 


0.17 


0.05 


0.00 


-0.96 


-0.01 


Sr4076 


1.235 


0.085 


-1.02 


-0.52 


0.50 


-0.05 


-0.39 


-0.06 


-0.90 


0.02 


-0.09 


-0.33 


-0.11 


Eu3970 


6.259 


0.080 


-0.67 


0.17 


-0.68 


-0.17 


-2.07 


0.85 


-1.45 


-0.16 


-0.01 


9.56 


-0.06 


Eu4592 


0.041 


0.083 


-1.05 


-0.12 


0.11 


-0.11 


-0.35 


-0.01 


0.36 


-0.02 


0.00 


0.50 


0.06 


on 


1.992 


0.189 


-0.54 


0.57 


0.19 


-0.02 


0.14 


-0.02 


0.71 


0.01 


0.00 


-0.11 


-0.01 


XX3580 


6.517 


0.282 


2.68 


4.89 


-5.50 


0.03 


-0.75 


0.06 


-0.12 


-0.29 


0.00 


-0.29 


0.15 


NH3360 


3.221 


0.285 


-0.05 


5.18 


-0.52 


-0.03 


-3.43 


-0.04 


-0.30 


-0.05 


0.00 


-0.20 


0.13 


NH3375 


4.688 


0.322 


-0.20 


4.51 


-0.87 


-0.03 


-0.71 


-0.04 


-0.04 


-0.05 


0.00 


-0.50 


-0.02 


Mg3334 


1.542 


0.105 


-0.43 


-0.22 


-1.14 


-0.04 


5.91 


-0.05 


0.14 


-0.04 


0.00 


-0.35 


-0.24 
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Table 1 Continued 



Index 


1 1 


V 




Mil 


re 




IN 1 




Zill 


or 


Tin 


Hill 


UpAz 


Co3701 


0.07 


0.39 


-0.43 


0.35 


-0.24 


0.52 


-0.36 


0.00 


0.00 


0.00 


0.00 


-0.05 


0.02 


Oo3840 


-0.79 


0.20 


-0.71 


0.30 


-0.87 


0.66 


-0.88 


0.00 


0.00 


0.00 


0.00 


0.00 


4.74 




-0.46 


0.22 


-0.28 


-0.05 


-0.82 


0.44 


-0.64 


0.00 


0.00 


0.00 


0.01 


-0.01 


-4.91 


Co7815 


-0.20 


0.00 


-0.01 


-0.02 


-0.33 


2.74 


0.04 


0.00 


0.00 


0.00 


0.00 


0.00 


-1.03 


Co8185 


-1.74 


0.05 


0.02 


0.01 


-0.15 


1.38 


-0.01 


0.00 


0.00 


0.00 


0.00 


0.00 


-4.53 


INlODDY 


-1.16 


0.10 


-0.04 


-0.14 


-6.07 


-0.64 


0.75 


0.00 


0.00 


0.00 


0.00 


0.01 


3.17 


JNioToU 


-2.63 


0.16 


-0.U7 


-U.ol 


-8.39 


f\ Art 
-U.4U 


1.7d 


0.00 


0.00 


u.uu 


0.00 


0.00 


2.31 


■\t; a or\o 

JNi4z9z 


2.47 


0.19 


0.25 


0.00 


-2.62 


0.07 


0.78 


-0.01 


0.00 


-0.01 


0.01 


0.00 


-1.81 


■\t; /( n 1 n 


-0.16 


-0.24 


-0.94 


-0.06 


-U.4z 


-0.03 


0.51 


0.00 


0.00 


-0.04 


0.03 


0.00 


0.61 


JNi497d 


1.07 


-0.01 


-0.34 


-0.02 


-3.22 


0.00 


0.62 


0.00 


0.00 


-0.06 


0.00 


0.00 


2.56 


"NT: cr cr no 

JNiooyz 


-1.18 


0.11 


0.03 


-0.05 


-0.28 


0.01 


0.59 


0.00 


0.00 


-0.01 


-0.03 


0.00 


-0.70 


Cu5217 


1.76 


0.04 


2.26 


-0. 14 
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APPENDIX B 
INDEX VELOCITY DISPERSION TRENDS 
The Virgo spectra (red symbols with error bars), the Virgo emission spectra (black cross 
symbols with error bars), the SDSS spectra (green symbols with error bars) and the Toloba 
spectra (black box symbols with error bars) are plotted. Single stellar population models 
(blue grid lines) from bottom to top of ages 17, 12, 8, 5, 3, 2, and 1.5 Gyrs and models 
from right to left of metallicities 0.5, 0.25, 0, —0.25, —0.5, —1, —1.5 and —2 (pink grid lines) 
are also plotted. 

The following 74 figures consist of one index plotted against o in the first panel, Mg h in 
the second and <Fe> in the third. In general the index is fitted with a single black best fit 
line for all data sets. However if one data set is significantly offset from the others then it 
will have its own fit line colored the same as that data set and the other two will be fit by a 
black line. A red line is the fit for the Virgo data both emission and non. A green line is 
the fit for the SDSS data. There are plots for which no trend could be discerned for those 
plots there is no line fit. 

The figures are arranged by order of atomic number starting with hydrogen; where there 
are more than one index for a given element or elements the indices are further arranged so 
that the bluest index is first. 

The results of this arrangement puts hydrogen first followed by C, N and O then the alpha 
elements followed by the iron peak elements. 

Note that in the Ha plot (Fig. 1) the Virgo emission galaxies are not plotted because they 
fall well below the range of the other measurements. 
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H beta vs Sigma 
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HdeltaA vs Sigma 
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CN 1 vs Sigma 



CN 1 vs Mgb 



CNl vs <Fe> 





0.3 




0.25 




0.2 




0.15 




0.1 




0.05 









-0.05 




-0.1 




Virgo Data 
Virgo Emission 
Sloan Data 
Toloba Data 



50 100 150 200 250 300 350 4QS 1 1.5 2 2.5 3 3.5 4 4.5 5 0.5 



Sigma in km s 



Mgb (A) 



1.5 2 2.5 3 
<Fe> (A) 



Fie. 9.- 



CN 2 vs Sigma 



CN2 vs Mgb 



CN2 vs <Fe> 



5 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 


-0.05 



Virgo Data 
Virgo Emission 
Sloan Data 
TolobaJ}ala 




0.35 

0.3^ 
0.25 


0.35/ 

yO,3 










0.2 


.y.^ 0.2 










0.15 


■A' o'i5 










0.1 


— 0.1 










0.05 



I ^i-'/ 0.05 










-0.05 


■ ■ ■ ' ■ / -0.05 





Sigma in km s 



50 100 150 200 250 300 350 4m 1 1.5 2 2.5 3 3.5 4 4.5 5 6.5 1 1.5 2 2.5 3 3.5 4 
--' Mg b (A) <Fe> (A) 



Fig. 10. 



59 



CN03862 vs Sigma 



CN03862 vs Mgb 



CN03862 vs <Fe> 



16 
14 
12 
10 
8 
6 
4 
2 




Virgo Data 
Virgo Emission 
Sloan Data 
Toloba Data- 



^16 
14 
12 
10 
8 
6 
4 
2 




50 100 150 200 250 300 350 4(0) 
-1 



1 1.5 2 



Sigma in km s 



16 
14 
12 
10 
8 
6 
4 
2 




2.5 3 3.5 4 4.5 5 0.5 

Mgb (A) 



1.5 2 2.5 
<Fe> (A) 



3.5 4 



Fiff. 11.- 



CN04175 vs Sigma 



CN04175 vs Mgb 



CN04175 vs <Fe> 




50 100 150 200 250 300 350 4QS 1 
Sigma in km s~^ 



1.5 2 2.5 3 3.5 4 4.5 5 

Mgb (A) 



6.5 



1.5 2 2.5 
<Fe> (A) 



Fie. 12.- 



C04685 vs Sigma 



C04685 vs Mgb 



C04685 vs <Fe> 



8 




Virgo Data 


8 




8 




7 


Virgo Emission 


7 












Sloan Data ! X 






■■ " 6 




6 


Toloba Data 


6 








5 








5 






. . , 5 






4 
3 








: 4 
3 






..l4 . 






2 








2 






' , u 2 






1 




i 1 


1 






i 

— 1 





























-1 








-1 






-1 







50 



100 150 200 250 300 350 m) 
-1 



Sigma in km s 



1 1.5 2 2.5 3 3.5 4 4.5 5 0.5 

Mgb (A) 



1.5 



2 2.5 
<Fe> (A) 



3.5 
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Fig. 20.- 
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Fig. 23.- 
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Fig. 29.- 



Sc4312 vs Sigma 
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APPENDIX C 
VIRGO CLUSTER GALAXY INDEX GRADIENTS 
The following 74 figures consist of index measurements of 18 Virgo cluster galaxies plotted 
as a function of log(r/ Re) where r is the radius and is the effective radius of each 
individual galaxy. In each figure there are 18 plots one for each galaxy and in each plot is a 
single black best fit trend line. 

The figures are arranged by order of atomic number starting with hydrogen, where there 
are more than one index for a given element or elements the indices are further arranged so 
that the bluest index is first. 

The results of this arrangement puts hydrogen first followed by C, N and O then the alpha 
elements followed by the iron peak elements. 

The five galaxies which show hydrogen emission in Ha are NGG 1400, NGG 2768, NGG 
3156, NGG 4278 and NGG 4486. 
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